
cR.-i/^s;2S3 

NASA CR - 165253 



NASA-CR- 165253 

OGI G&QCI 


ION AND ADVANCED ELECTRIC THRUSTER RESEARCH 


PREPARED FOR 
LEWIS RESEARCH CENTER 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


GRANT NGR-06-002-112 


I-. wJii' I 

MAR 9 1981 

U.(, -;r.!j!:ARCM CENTSfi, 
LIBRARY, NASA' 

.'wwifi-tMN, 'WiiwtenAi 

Annual Report 
December 1 980 
Paul J. Wilbur 

Department of Mechanical Engineering 
Colorado State University 
Fort Collins, Colorado 



1. Report No. 2. Government Accession No. 

NASA CR 165253 

3. Recipient's Catalog No. 

4. Title and Subtitle 

Ion and Advanced Electric Thruster Research 

5. Report Date 
Dec. 1980 

6. Performing Organization Code 

7. Author(s) 

Paul J. Wilbur 

8. Performing Organization Report No. 

10. Work Unit No. 

g. Performing Organization Name and Address 

Department of Mechanical Engineering 
Colorado State University 
Fort Collins, Colorado 80523 

11. Contract or Grant No. 
NGR-06-002-112 

13. Type of Report and Period Covered 

Contractor Report 
Dec. 1 . 1979 - Dec. 1 , 1980 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 



14. Sponsoring Agency Code 

1 15. Supplementary Notes 


Grant Monitor - William Kerslake, NASA Lewis Research Center, Cleveland, Ohio 44135. 


A phenomenological model of the orificed, hollow cathode based on the field-enhanced, 
thermionic mechanism of electron emission is presented. Experimental results obtained on a 
mercury hollow cathode are presented and they are seen to support the validity of the model . 

High frequency oscillations associated with the orificed, hollow cathode are shown to be a 
consequence of current flow through the cathode orifice, A procedure for Langmuir probing 
of the hollow cathode discharge and analyzing the resulting probe characteristics is discussed. 
The results of sputter yield measurements made for molybdenum, tantalum, type 304 stainless 
steel and copper surfaces being bombarded by low energy argon or mercury ions are presented. 

The experimentally determined effects of nitrogen and alternated copper layers on the sputter 
yields of molybdenum, tantalum and 304 stainless steel are also discussed. A dynamic model of 
electrothermal rocket and ramjet thrusters is developed. The gross performance of these 
deyjce^ js compared to that of an electromagnetic gun for the case of a high acceleration, earth 
launch mission. The theoretical performance of electrothermal rockets and ramjets is shown to 
be comparable to that of the electromagnetic gun. The experimental apparatus that will be used 
to measure the divergence and deflection characteristics of ion beamlets as a function of grid 
aperture misalignment is discussed. Preliminary results obtained with this apparatus are 
presented. The experimentally determined operating characteristics of two hollow cathodes 
connected in parallel to supply electrons to a thruster discharge chamber are presented. 

Results show these cathodes can be started and forced to share equally in the required electron 
emission. A model describing the effect of primary electron energy on the ion velocity needed 
for a stable sheath (Bohm velocity) is presented. 


17. Key Words (Suggested by Author(s)) 

Electrostatic Thruster 
Hollow Cathode 
Electrothermal Thruster 

18. Distribution Statement 

Unclassified - Unlimited 

19. Security Classif. (of this report) 
Unclassified 

20. Security Classif. (of this page) 
Unclassified 

21. No. of Pages 
130 

22. Price* 


* For sale by the National Technical Information Service, Springfield, Virginia 22161 


NASA-C-168 (Rev. 10-75) 


i 






TABLE OF CONTENTS 


Topic Page 


Abstract i 

Ion Extraction from a Plasma 1 

Baffle Aperture Design Study of Hollow Cathode 

Equipped Ion Thrusters 1 

A Phenomenological Model Describing Or if iced. 

Hollow Cathode Operation 3 

Background 3 

A Phenomenological Model 6 

Comparison with Experiment 10 

Plasma Density in the Ion Production Region 13 

Surface Work Function and the Emission Mechanism 23 

Future Work 28 

Conclusions 30 

Hollow Cathode Plasma Oscillations . 32 

Experimental Results 32 

Discussion 37 

Conclusions 39 

Probe Measurements inside the Orificed, Hollow Cathode 40 

Spherical Probe Design 40 

Data Analysis 42 

Probe Contamination 46 

Plasma Oscillations 49 

Additional Considerations 50 

Conclusions 52 

Sputtering Rates of Discharge Chamber Materials 53 

Apparatus and Procedure 54 

Results 59 

SERT II Thruster Tests 68 

Electrothermal Rocket 74 

Constraints at End of Launch Phase 77 

Results 81 

Conclusions 88 


ii 



Table of Contents (Continued) 

Ion Beamlet Vectoring by Grid Translation 89 

Apparatus and Procedure 90 

Results 94 

Conclusion 97 

Multiple Cathode Study 98 

Introduction 98 

Apparatus ..... 98 

Discussion 100 

Concluding Remarks 105 

Modified Bohm Velocity for Plasmas Containing 

Primary Electrons and Doubly Charged Ions, . . , . 107 

References 113 

Appendix A - Langmuir Probe Sweep Circuitry 116 

Appendix B - Electrothermal Launch System Dynamics 118 

Distribution List 124 



LIST OF FIGURES 


Figure No. Title 

1. Typical Orificed, Hollow Cathode. . 

2. Test Cathode for Current Accounting 

3. Typical Current-Voltage Characteristic 

4. Effect of Internal Cathode Pressure on Current 

Distribution 

5. Effect of Temperature on Insert Work Function 

6. Insert with Test Patch 

7. Apparatus for In-situ Measurement of Insert Work 

Function 

8. Quartz Test Cathode 

9. Probe Floating Potential , Fluctuations 

(I|^ = 3.5A, = 16.6v, m = 100mA, Spot Mode) 

10. Probe Floating Potential , FI uctuations 

(I^ = 2.3A, = 21 .8v, m = 100mA, Plume Mode) 

11. Probe Floating Potential ^Fl actuations 

(l^ = 5.3A, = 17.5v, ifi = 100mA, Spot Mode) 

12. Probe Floating Potential -El ectron Collection 

at Orifice Plate (Ir = 2.3A, V. = 12.lv, m = 100mA, 

Plume Mode) 

13. Spherical Probe Design 

14. Typical Probe Trace 

15. Current Density Profile 

16. Sputter Erosion Test Apparatus 

17. Argon Ion Sputtering Yields 

18. Mercury Ion Sputtering Yields 

19. Effect of Nitrogen on Argon Sputtering Yields 

20. Effect of Nitrogen on Mercury Ion Sputtering Yields . . . 

21. SERT II Neutralizer Keeper Current Collection 

Characteristic - No Neutralizer Flow 

iv 


Page 

4 

12 

16 

20 

25 

26 

29 

33 

35 

35 

36 

36 

41 

43 

55 

56 
60 
61 

63 

64 

70 



22. SERT II Beam Plasma Properties 71 

23. SERT II Test Circuit Diagram 73 

24. Electrothermal Thruster Schematics 76 

25. Launch Phase Nomenclature 78 

26. Velocity Loss Through Atmospheric Drag 80 

27. Velocity Profile Comparison 82 

28. Acceleration Profile Comparison 84 

29. Mass Profile Comparison 86 

30. Power Profile Comparison 87 

31. Deflecting Beamlet Nomenclature 91 

32. Deflecting Beamlet Test Apparatus 92 

33. Beamlet Deflection Characteristic 95 

34. Beamlet Divergence Characteristic 96 

35. 15 cm Multi pole. Multicathode Thruster 99 

36. Multicathode Thruster Circuit Diagram 101 

37. Location of Tuning Resistors for Control of 

Relative Cathode Emission 103 

38. Cathode Emission Control Using Tuning Resistors 104 

39. Control Parameter Profiles 106 

A1 . Probe Voltage Sweep Circuit 117 


V 



LIST OF TABLES 


Table No. Title Page 

I. Cathode Current Accounting 14 

II. Calculated Plasma Densities 14 

III. Effect of Insert Position on Cathode Densities 18 

IV. Effect of Point of Electron Collection on Cathode 

Operating Parameters 18 

V. Sputter Erosion Characteristics of Multilayered 

Test Samples in an Argon Ion Beam 65 

VI. Sputter Erosion Characteristics of Multilayered 

Test Samples in a Mercury Ion Beam 67 

VII. SERT II Thruster Flow Rate Data 

(High Voltage Off) 68 


vi 



ION EXTRACTION FROM A PLASMA 


Graeme Aston 

An experimental investigation of the physical processes governing ion 
extraction from a plasma was completed during this grant period. The re- 
sults of this work are described in Reference 1, so they will not be 
reproduced here. In Reference 1 the screen hole plasma sheath of a multi- 
aperture ion accelerator system is defined by equipotential plots for a 
variety of accelerator system geometries and operating conditions. A 
preacceleration region, which might be called the sheath, having a thickness 
of at least fifteen Debye lengths is shown to be typical. The electron 
density variation within the sheath is shown to satisfy a Maxwel 1 -Bol tzmann 
density distribution at an effective electron temperature dependent on the 
discharge plasma primary-to-Maxwell ian electron density ratio. Plasma ion 
flow up to and through the sheath is shown to be predominately one-dimensional 
with the ions entering the sheath at a modified Bohm velocity. Low values 
of the screen grid thickness-to-screen hole diameter ratio are shown to 
give good ion focusing and high extracted ion currents because of the effect 
of screen webbing on ion focusing. 

BAFFLE APERTURE DESIGN STUDY OF HOLLOW CATHODE EQUIPPED ION THRUSTERS 

John Brophy 

A simpl e, theoretical model which can be used as an aid in the design 
of the baffle aperture region of a hollow cathode equipped ion thruster 
was developed during this grant period. The results of this study are 
described in detail in Reference 2 and its contents will also not be 
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reproduced here. The document contains an analysis of the ion and electron 
currents in both the main and cathode discharge chambers. From this 
analysis a model of current flow through the aperture, which is required 
as an input to the design model, is developed. This model is verified 
experimentally. The dominant force driving electrons through the aperture 
is shown to be the force due to the electrical potential gradient. The 
diffusion process is modeled according to the Bohm diffusion theory. A 
number of simplifications are made to limit the amount of detailed plasma 
information required as input to the model and to facilitate the use of the 
model in thruster design. This simplified model can be used to predict the 
magnetic field intensity required to produce a given current through the 
baffle aperture in terms of the anode/keeper voltage difference and the 
baffle aperture/magnetic field geometric parameters. The magnetic field 
intensities computed from the model are shown to be remarkably consistent 
with experimentally measured values obtained with a given thruster geometry 
over substantial changes in operating conditions. Dramatic changes in the 
thruster geometry can introduce effects that cause a coefficient in the 
equation describing the model to vary by as much as a factor of two. The 
design usefulness of the model may be limited by this factor of two un- 
certainty, which exists when the thruster geometry is altered substantially, 
and by the accuracy to which the plasma parameters required as inputs to 
the model can be specified. 
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A PHENOMENOLOGICAL MODEL DESCRIBING 
ORIFICED, HOLLOW CATHODE OPERATION 

Dan Siegfried 

Because of the importance of the orificed, hollow cathode to ion 
thruster design and performance, there has been considerable interest in 
recent years in understanding hollow cathode physical processes and in 
developing a model for cathode operation. Numerous experiments have been 
performed under the support of this grant which have provided important 
physical data regarding hollow cathode operation. These data have in- 
cluded plasma properties within the cathode, insert temperature profiles, 

internal cathode pressure, and emission current density profiles over a 

3 4 

wide range of cathode operating conditions. ’ The physical understanding 
gained through this work has led to the development of a simple phenomeno- 
logical model to describe cathode operation. This report will discuss our 
current understanding of the physical processes which take place inside 
the hollow cathode and describe the proposed model. Work on the model has 
suggested additional experiments for verifying the model and these will 
also be discussed. 

Background 

The orificed, hollow cathode in its simplest form is shown in Fig. 1. 
The cathode consists of an outer refractory metal tube, usually tantalum, 
covered on its downstream end by an orifice plate usually made of thoriated 
tungsten. The cathode also normally incorporates a refractory metal insert 
either coated or impregnated with such chemicals as various barium com- 
pounds, which aid the emission process by reducing the work function of 




Fig. 1. Typical Orificed, Hollow Cathode 
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the insert surface. Cathodes used in ion thrusters typically have tube 
inside diameters of a few millimeters and orifice diameters of a few 
tenths to one millimeter. The insert length is usually a few tube dia- 
meters. The electron current is collected by an anode biased positive of 
the cathode. Thruster applications also utilize a small secondary anode, 
called a keeper, which acts to maintain a stable discharge in the event 
of fluctuations in discharge chamber conditions. Mercury vapor flows 
through the cathode in the direction shown in the figure. An external 
heater (not shown) is normally used to heat the cathode as an aid to 
starting the emission process. However, the process is self-sustaining 
once the cathode discharge is established and the heater can then be turned 
off. 

The earlier experimental studies mentioned in the introduction have 
led to a number of important conclusions which have helped to define the 
nature of the hollow cathode operation. These will be summarized very 
briefly here because they form the basis of the proposed phenomenological 
model. Langmuir probe measurements made within the cathode indicated that 

3 

there is a very intense plasma discharge set up within the cathode. For 

emission currents of a few amperes, typical plasma conditions found in the 

vicinity of the orifice plate were a plasma density of a few times 
14 -3 

10 cm , a plasma potential of 8 volts, and an electron temperature of 
0.8 eV. A series of experiments using a quartz body test cathode to 
isolate various cathode surfaces electrically was used to pinpoint the 
emission site in the cathode and to determine emission current profiles for 

4 

the insert. That study showed that the insert accounts for approximately 
85% of the total emission current and that the emission region is localized 
on the last few millimeters at the downstream end of the insert. This same 
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series of tests also provided insert temperature profiles over a wide 
range of operating conditions. As expected the highest insert temperature 
was found to correspond to the location of the emission region, with 
temperatures on the order of 1000°C being typical for an emission current 
of a few amperes. Furthermore, it was found that the orifice diameter 
and mass flow rate were important only to the extent that they affected 
the internal cathode pressure, and that the cathode pressure was an im- 
portant parameter affecting the emission temperature. The product of 
cathode internal pressure and emission surface diameter was suggested as 
an important parameter in determining the emission location and keeper 
voltage. This parameter is proportional to the ratio of the insert diameter 
to the electron mean free path, and therefore, is also proportional to the 
probability of an inelastic collision of an electron as it makes a single 
traverse of the cathode diameter. The simultaneous measurement of emission 
current and emission surface temperature in these studies provided suffic- 
ient information to facilitate estimation of the emissive surface work 
function. This provided reasonable evidence that the emission mechanism of 
the cathode could be that of field-enhanced, thermionic emission. 

A Phenomenological Model 

The following physical description of orif iced, hollow cathode operation 
has emerged from these studies. The cathode orifice maintains a high 
neutral density inside the cathode by restricting the propellant flow and 
it also provides a current path to the downstream discharge. The electron 
emission comes uniformly from a narrow {'^ 2 mm) band on the downstream end 
of the insert. The electrons are produced at the surface of the insert by 
field-enhanced, thermionic emission (the very strong electric field is a 
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consequence of a very dense plasma and the resulting potential drop across 
a very thin plasma sheath). The electrons produced at the insert surface 
are accelerated across the plasma sheath by a potential of 8 to 10 V. 

Since the mean free path for inelastic collisions of these energetic 
electrons is on the order of the internal cathode diameter, the "ion 
production" region can be idealized to be the volume circumscribed by 
the emitting region of the insert. The dense internal plasma is estab- 
lished by the ionization taking place in this region. Ions produced in 
this volume diffuse out of it at the Bohm velocity. These ions strike 
the insert surface with sufficient energy to heat it to the temperature 
necessary to provide the required electron emission. The emission sur- 
face temperature, however, is determined not only by the emission current 
but also by the local plasma properties. 

The plasma properties in the ion production region are coupled into 
the problem by the energy balance at the insert surface in the following 
manner. The plasma properties determine the ion flux and therefore the 
energy input to the emission surface. For a given emission current, the 
surface temperature is determined by the energy balance which demands that 
the thermal losses from the surface due to electron production, radiation 
and conduction are balanced by the energy input from the ion flux. The 
plasma properties also affect the required emission temperature because 
they determine the magnitude of the electric field at the emission surface 
and thereby the degree of field-enhancement in the emission process. 
Therefore, for a given emission current the surface temperature and plasma 
properties must be consistent to the extent that they satisfy the energy 
balance at the surface. All cathode surfaces which contact the plasma 
receive ion currents proportional to the Bohm velocity and the plasma 



8 


density adjacent to the surface. Electron emission, on the other hand, 
can be assumed to come only from the 2 mm band on the downstream end of 
the insert. The total emission current* from the cathode is equal to the 
sum of the ion currents to the various cathode surfaces and the current 
of the emitted electrons. 

Certain aspects of this phenomenological model can be expressed 
analytically in a simple form which will allow comparison with experimental 
results. The plasma density adjacent to a particular surface (n) can be 
calculated based on the Bohm condition using 


®'^Bohm 

I. 

1 

■"e‘ 

h 

(1) 

Ae 

■V- 


i 


where is the ion current to the surface, A is the surface area, is 

the Maxwellian electron temperature (°K), m^. is the ion mass, e is the 

electronic charge, and k is Boltzmann's constant. For an electron emitting i 

surface the measured current to the surface is determined by both collected 

ions and emitted electrons, so that the total current density to the ' 

surface is 

Hotal " '^i ^ ^e 

where j^. is the ion current density and j^ is the electron emission current 
density. The ratio of ion to electron currents can be estimated from an 
energy balance on the emitting surface. In such a balance the ion heating 

power is equated to power conducted and radiated from the surface plus the j 

5 ' 

power required to boil off electrons. The equation describing this is 

* Total emission current (l^) as used here refers to the total power 
supply current, that is the keeper current and the anode current. 
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J'e *e ^ (''c * ''1 - 


(3) 


where is the effective work function of the surface, q is the thermal 
heat flux away from the surface, is the potential drop across the plasma 
sheath, V^. is the ionization potential, and c|)^ is the work function of the 
surface material (a material property). Equations 2 and 3 can be combined 
to give the electron emission current density from the surface. 


'total 


- a q 


1 + a 


(4) 


where 

a MVc V, - . 

In general the thermal loss is a function of the surface temperature and 
the cathode thermal design. For the quartz test cathode which will be 
discussed later most of the thermal loss is due to radiation from the outer 
surface of the insert to father cold external surfaces, and can be estimated 
from 

q ~ eaT** (5) 

where e is the emissivity 0.5 for tantalum), a is the Stefan-Bol tzmann 

constant, and T is the surface temperature. Emission from the surface is 

5 

assumed to be given by the Schottky equation for field-enhanced, thermionic 
emission _ 

Je ' V [- 

where =120 A/cm^K^ and the other parameters are as previously defined. 
The average effective work function is given by 
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where is the permittivity of free space and E, the electric field 
adjacent to the surface, can be estimated using 
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Here the factor of 4/3 comes from Child's Law considerations and the 
sheath thickness is estimated as one Debye length (Xg). 

This model as it stands does not contain sufficient detail to make 
a priori predictions of plasma properties or their spatial variation. The 
model, therefore, cannot explain or predict the location at which the 
emission will take place or the performance, i.e. power requirements, of 
the cathode. This would require a more detailed modeling of the plasma 
processes in the ion production region. However, the model can provide an 
estimate of the insert temperature and this is a critical parameter in 
determining both the cathode lifetime and performance. This can be done 
for example, by picking an electron temperature and the plasma potential. 
These properties have been measured experimentally and found to be rather 
insensitive to operating conditions with typical values of '\'0.8 eV and 'V/8.0 
volts respectively for a cathode operating at a few amperes of discharge 
current. Using these values together with a specified surface work function 
((j)^) and the desired emission current (j-^otal^’ ^ through 8 can be 
solved to provide the emission surface temperature (T). Note that the 
solution scheme would generally be an iterative one requiring an initial 
guess of a value for either j^ or T. 

Comparison with Experiment 

The model has been found to be consistent with the experimental results 
obtained to date. As a result of the work on the model a number of addi- 
tional experiments were proposed. One such experiment was an overall 
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current balance on all of the cathode surfaces which can contribute to 
the total emission current. Such an experiment would provide a means of 
testing the assumptions related to the ion production region such as the 
extent of the region, the applicability of the Bohm velocity, and the 
energy balance at the insert surface. 

The test cathode designed and built to facilitate the measurement of 
the currents necessary to perform the overall current balance is shown in 
Fig. 2. The cathode was fabricated from a quartz tube such that all of 
the internal surfaces were covered with tantalum foil and each surface of 
interest was isolated as a separate electrode. The insert was made 2.2 mm 
long to test the assumption regarding the location and extent of the emis- 
sion region. Both the insert and the cathode tube were coated with the low 
work function chemical R-500*. This cathode had an inner diameter of 3.9 mm 
and an orifice diameter of 0.76 mm. The cathode discharge was coupled to a 
cylindrical anode which was completely enclosed within a stainless steel 
ground tube covered on its downstream end by a fine mesh screen. The 
current from each of the surfaces shown in Figure 2 was measured separately 
for total emission currents of 1.3, 2.3, 3.3 and 4.3 A at a flow rate (m) 
of '\'100 mA. Data were also collected for an emission current of 3.3 A at 
internal cathode pressures (P) ranging from 1.3 to 5.8 torr. Internal 
pressures were measured using a U-tube manometer as described in the report 

4 

of an earlier experimental study. The internal pressure was adjusted by 
venting a portion of the mercury flow into the bell jar using a special 
throttle valve also described in that report. The results of these 

* A double carbonate mixture - (Ba/Sr) CO^ - manufactured by the J. T. 

Baker, Co.; Phil 1 ipsburg. New Jersey. 
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Fig. 2. Test Cathode for Current Accounting 
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experiments are shown in Table I where the currents - Ii, l 2 » etc. - 
are those from the surfaces identified in Fig. 2 and suggested by the 
diagram attached to Table I. The current is the sum of these numbered 
currents. It appears that all of the currents necessary to perform the 
overall cathode current balance are accounted for because the total of 
the currents (1^) agrees with the total emission current (I^ " 
within a few percent. 

Plasma Density in the Ion Production Region 

One of the purposes of this test was to determine the fraction of the 
total emission current resulting from the ion current returning to cathode 
potential surfaces. The insert segment, which in this test was 2.2 mm long, 
accounted for '\^85% of the total emission current except at pressures below 
about 3 torr. According to the simple model proposed earlier, this insert 
current is equal to the sum of the electron emission current (given by the 
Schottky equation - Eq.6) and the ion collection current (based on the 
Bohm velocity - Eq.1). The ratio of the ion current to the electron current 
for the insert was estimated from the insert energy balance (Eq. 4) and was 
found to be in the range of 0.17 to 0.24 for the operating conditions of 
Table I. Using this ratio together with the measured insert current (I 2 ) 
one can estimate the ion current to the insert. 

The cathode model which has been proposed, also assumed that electron 
emission took place only at the insert surface and that the current contri- 
bution of the other cathode surfaces is due solely to the ion currents to 
these surfaces. This assumption was tested experimentally by biasing each 
of the electrode surfaces (except the insert) with respect to ground and 
observing the effect on the collected current. A typical voltage current 



TABLE I 


m(mA) 

'^(torr) 

VA) 

Va(V) 

I,(A) 

Cathode 

Vk(V) 

Current Accounting 

I^(A) I^CA) 13 (A) 

(Cathode ( Insert) (Upstm. 
tube) Orifice 

Plate) 

102 

3.6 

1 .0 

15.8 

0.3 

13.9 

0.03 

1.10 

0.10 

102 

4.6 

2.0 

14.2 

0.3 

12.8 

0.05 

1 .99 

0.16 

106 

5.5 

3.0 

15.2 

0.3 

14.0 

0.07 

2.82 

0.23 

105 

6.2 

4.0 

15.7 

0.3 

14.2 

0.09 

3.63 

0.29 

no 

5.8 

3.0 

15.1 

0.3 

15.0 

0.05 

2.82 

0.24 

78 

4.1 

3.0 

15.1 

0.3 

15.5 

0.14 

2.78 

0.24 

62 

3.2 

3.0 

15.2 

0.3 

16.2 

0.25 

2.72 

0.23 

35 

1.8 

3.0 

19.0 

0.3 

17.4 

1.25 

1.76 

0.18 

25 

1 .3 

3.0 

24.1 

0.3 

17.1 

1.90 

1 .23 

0.14 




TABLE II 



Calculated 

Plasma Densities 

Cathode 

Internal 

Pressure 

(torn) 

Total Emission 

Current 

Ie(a) 

Insert 

Density 

n 2 (cm-^) 

Density 
Immediately 
Upstream of 
Orifice Plate 
n3(cra-3) 

3.6 

1.3 

0.79 X lOl** 

0.91 X 101“ 

4.6 

2.3 

1.21 X lOl"* 

1 .46 X 101“ 

5.5 

3.3 

1.58 X 101“ 

2.10 X 101“ 

6.2 

4.3 

1.96 X 101“ 

2.64 X 10'“ 

5.8 

3.3 

1.58 X 101“ 

2.19 X 101“ 

4.1 

3.3 

1.58 X 101“ 

2.19 X 101“ 

3.2 

3.3 

1.55 X 101“ 

2.10 X 101“ 

1.8 

3.3 

1.21 X 101“ 

1.64 X 101“ 

1.3 

3.3 

0.98 X 101“ 

1 .28 X 101“ 


).r“ a" 1 




I,(A) Is(A) I, (A) Ij(A) 
(Orifice)(Downstm. (Ground (*C) 

Orifice Screen) 


Plate 

0.04 

3.0 

0.002 

1 .27 

965 

0.09 

0.01 

0.004 

2.30 

1009 

0.17 

0.02 

0.007 

3.31 

1039 

0.26 

0.025 

0.013 

4.30 

1070 

0.17 

0.015 

0.007 

3.30 


0.11 

0.01 

0.007 

3.29 

— 

0.09 

0.01 

0.006 

3.31 


0.05 

0.01 

0.006 

3.26 

— 

0.04 

-vO 

0.006 

3.32 

-- 


Density 

Within 

Orifice 

n4(cm-i) 

Density 

Immediately 

Downstream 

of Orifice 

Plate 

H5(cm"i) 

0.75 X 101“ 

^*0 

1 .70 X 101“ 

2.39 X IOI 2 

3.21 X 101“ 

4.78 X IOI 2 

4.90 X 101“ 

5.98 X I 012 

3.21 X 101“ 

3.58 X IOI 2 

2.08 X 101“ 

2.39 X I 012 

1.70 X 101“ 

2.39 X IO 12 

0.94 X 101“ 

2.39 X IO 12 

0.75 X 101“ 

^»o 
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trace obtained from these tests is shown in Fig. 3. Since the curve is 
approximately horizontal where it crosses the y-axis (i.e., at cathode 
potential), the current does appear to be primarily a collected ion current. 
This test demonstrates that the surfaces collect negligible electron current 
but does not prove whether or not electrons are being emitted from the 
surface. It is believed, however, that the surface work functions are 
sufficiently high to preclude significant electron emission from these 
surfaces at their prevailing temperatures. Electron production due to 
photoemission or from ion-induced, secondary emission have also been 
estimated analytically and found to represent insignificant components of 
the measured current. On the basis of these considerations currents I^, 

I3, 14 , I5 and Ig-can be assumed to give the ion current to each corresponding 
cathode surface. Adding these to the insert ion current obtained in the 
preceding paragraph one obtains the total ion current to cathode surfaces. 
The result of this is that in general about 32% of the total emission cur- 
rent (I^) is due to ions returning to cathode potential surfaces. About 
17% of this total emission current is ion current to the insert and the 
remainder (15%) to the other cathode surfaces. This distribution of currents 
appears to be valid over the full range of emission current investigated. 

As the cathode pressure is reduced the cathode tube begins to emit(Ii in 
Column 7 of Table I), but the basic distribution of ion and electron currents 
suggested by this analysis is still considered to be valid (i.e., '^M% of 
the emission current is actually ion current to the emitting region and 
^^15% is ion current to other cathode surfaces). 

* The ion current to the insert is~20% of the total insert current and the 
total insert current is 85% of the total emission current. 
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Fig. 3. Typical Current-Voltage Characteristic 


Knowing the ion current to each surface makes it possible to estimate 
the plasma density adjacent to each surface. This is done by using Eg. 1 , 
assuming that the surface is in a region of uniform plasma density and 
that the ions leave the plasma at the Bohm velocity based on an electron 

3 

temperature obtained in previous experiments of 'vO.8 eV. The calculated 

plasma densities adjacent to each surface are shown in Table II. The most 

notable result shown in this table is that the plasma densities adjacent 

to the insert (n2) and upstream side of the orifice plate (03) are close 

to each other over the full range of experimental conditions. This suggests 

that the ionization is, in fact, taking place in a narrow slice bounded by 

the insert segment and the upstream side of the orifice plate. It also 

suggests that the plasma density is relatively uniform in that region, and 

that the energy balance used to estimate the ion current to the insert 

yields self-consistent results. It should also be noted that the plasma 

densities shown in Table II for the region adjacent to the insert (n2) are 

in reasonable agreement with those found from Langmuir probe measurements 

3 

made in earlier studies under similar conditions. 

An experiment was also conducted using the cathode configuration of 
Fig. 2 but with the insert moved back '\^2 mm from the orifice plate. The 
plasma densities adjacent to the various cathode surfaces are shown in 
Table III for both the case where the insert is moved back from the orifice 
plate and the case where it is in the normal position. A comparison of 
the results shows that while the plasma density is roughly the same 
adjacent to both the insert and the orifice plate with the insert in the 
normal position, this is definitely not the case when the insert is moved 
upstream. With the insert upstream, the density adjacent to the orifice 
plate is only about one-third of that in the insert region but it jumps 
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Effect of Insert Position on Cathode Densities 


Insert 

Cathode 

Total 

Insert 

Density 

Density 

Density 

Position 

Internal 

Emission 

Density 

Immediately 

within 

Immediately 


Pressure 

Current 

ni(cm"3) 

Upstream of 

Orifice 

Downstream of 


(torr) 

Ip(A) 

Orifice Plate 

n3(cm-3) 

Orifice Plate 


h 


n2(cm"^) 

n4(cm"^) 


Normal 

5.8 

3.3 

1 .6 X 101“* 

2.2 X 101“+ 

3.2 X lOl** 

3.6 X 1012 

Upstrm 
2 mm 

5.8 

3.3 

1 .2 X 101'+ 

0.36 X 101'+ 

3.9 X 101*+ 

4.3 X 1012 


Anode 

Surface 

Cathode 

Internal 

Pressure 

(torr) 

Effect of Point 

Anode 

Current 

Ia(A) 

TABLE IV 

of Electron Collection on Cathode Operating Parameters 
(Insert 2 mm Upstream of Orifice Plate) 

Anode Keeper Keeper Insert 

Voltage Current Voltage Density 

V^(V) I|^(A) (V|^(V)) ni(cm-3) 

Density 
Immediately 
Upstream of 
Orifice Plate 
n2(cm“3) 

Normal 

4.6 

■ 3.0 

16.6 

0 

6.0 

1 .1 X 101*+ 

3.1 

X 1013 

Downstream 
Side of 
Orifice Plate 

4.5 

3.0 

12.3 

0 

5.5 

1 .1 X 101'+ 

3.1 

X 1013 

Upstream 
Side of 
Orifice Plate 

2.6 

3.0 

9.8 

0 


X 

o 



Normal 

2.6 

3.0 

17.8 

0 

6.3 

1 .1 X 101'+ 

3.1 

X 1013 


19 


by an order of magnitude as one proceeds into the orifice. This provides 
additional support for the assumption that the ion production region is 
that region circumscribed by the emitting portion of the insert and follows 
the insert when it is moved upstream. In addition, this suggests that the 
orifice is also an ion production region since with the insert moved upstream 
the density decreases in the region between the insert and the orifice plate 
presumably due to wall losses and then increases again within the orifice. 

The effect of internal cathode pressure on the current distribution of 

the cathode with the insert in its normal position is shown in Fig. 4. In 

this figure the measured current to the internal surfaces (from Table II) 

including the orifice bore is plotted as a function of internal cathode 

pressure for the cathode operating at a constant total emission current of 

3.3 A. This figure indicates that at pressures above 'v3 torr essentially 

all of the electron emission comes from the insert. Recall that according 

to the proposed model the electron emission current density (j^) from an 

emitting surface is proportional to the measured current density (j^otal^ 

or in this case I 2 /A 2 . At pressures below '^^3 torr, the insert emission 

current begins to drop off and the current from the cathode tube (Curve Ij 

in Fig. 4) begins to increase. This suggests that the electron emission 

region, and therefore the ion production region, spreads upstream along the 

cathode tube as the pressure is decreased below '^'3 torr. This increase in 

emission area in the upstream direction as pressure is decreased has been 

4 

seen in earlier experiments and previously reported. The curves in Fig. 4 
for measured current to the orifice plate (1 3 ) and the orifice bore (I 4 ) 
are also interesting in that they display different sensitivities to internal 
pressure. Assuming that there is no electron emission from these surfaces, 
it was argued earlier that these currents must be purely ion currents and 



MEASURED ELECTRODE CURRENT (A) 



Fig. 4. Effect of Internal Cathode Pressure on Current Distribution 
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as such are proportional to the plasma density adjacent to the respective 
surfaces (values of these densities are indicated in columns 4 and 5 of 
Table II). The plasma density adjacent to the orifice plate, as indicated 
by Curve I 3 (Fig. 4) follows the same behavior as the insert Curve I 2 . It 
is approximately constant at pressures above '^3 torr and decreases for 
pressures below that. The interesting thing, however, is that the ion 
current to the orifice bore (Curve I4, Fig. 4), and therefore the plasma 
density in that region (Column 5 of Table II), are linear functions of 
pressure. This is a different pressure dependance than that of plasma 
densities adjacent to the other internal surfaces. This suggests that the 
orifice plasma is isolated from the internal cathode plasma by either the 
gas dynamic effect of acceleration through the orifice or by the double- 
sheath commonly associated with a discharge constriction such as the orifice. 

The fact that the orifice plasma appears to be isolated from the in- 
ternal cathode plasma by some sort of a boundary or sheath was also 
illustrated by another experiment. Since the surfaces of the cathode were 
electrically isolated from one another, it was possible to bias each of 

them sufficiently positive with respect to the insert so they would collect 

all or a portion of the emission current. Such an experiment was conducted 
with rather interesting results. The cathode was operated at a total 
emission current of 3.0 A, without any keeper current and with a mass flow 
rate of 'v^lOO mA. The insert just happened to be moved 2 mm upstream of 

the orifice plate in this test but this is not felt to be significant. The 

results of collecting the emission current on the downstream side of the 
orifice plate (Surface 5 on Fig. 2) and on the upstream side of the orifice 
plate (Surface 3) are shown in rows two and three respectively of Table IV. 
Row one of this table gives the characteristics of the normal operating mode 
where the electrons are collected on an anode downstream of the keeper. The 
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most notable comparisons that can be seen in Table IV are between the 
potentials of the surface collecting the current (Column 4) and the 
internal pressures within the cathode (Column 2). Collecting the current 
at the downstream side of the orifice plate drops the anode voltage from 
a normal value of 16.6 v to 12.3 v but has only a small effect on the 
internal pressure. Collecting the emission current at the upstream side of 
the orifice plate, however, not only drops the anode voltage to 9.8 volts 
(this should be near the internal plasma potential) but also causes the 
internal cathode pressure to drop by nearly a factor of two. 

The fourth row in Table IV gives the results for the cathode coupled 
to the normal downstream anode but with a reduced flow so that the pressure 
is the same as in the case for results in the third row (collection at the 
upstream side of the orifice plate). Note that the calculated plasma 
densities adjacent to the insert are not significantly affected by any of 
these changes. The above results are consistent with earlier experimental 
results which have suggested that the orifice channel is a region that 
sustains strong potential gradients possibly in the form of a double sheath 
at the upstream end. These potential gradients could be responsible for 
the ion production which results in the high plasma densities in this region 
and are most certainly an indication of the plasma heating taking place 
there. This heating and the fact that ions in this region see an electric 
field driving them upstream may account for the reduction in upstream 
stagnation pressure which was observed when the heating in the orifice was 
removed by collecting the emission current at the upstream side of the 
orifice plate. It is also suggested that since the orifice provides the 
conduction path between the internal and external plasmas, that the orifice 
region relies on the potential drop across it to provide the energy for ion 



23 


production in the orifice region. When the emission current is collected 
downstream of the orifice, the resulting increased plasma density in the 
orifice region is needed to maintain current continuity in the presence of 
the large change in cross sectional area of the conduction path. 


Surface Work Function and the Emission Mechanism 

The most basic assumption of the proposed model is that the emission 
mechanism is that of field-enhanced, thermionic emission (described by 
Eqs. 6 and 7). The ability to verify this assumption will probably be the 
major test of the validity of the model. So far the experimental evidence 
points to field-enhanced, thermionic emission as being the dominant surface 
phenomenon producing electrons. This evidence was based at least in part on 
the ability of the Schottky equation (Eqs. 6 and 7) to predict reasonable 
values for estimates of the surface work function (fg).^ Since the model 
through Eq. 4 provides a means for determing the electron emission portion of 


the measured insert current, that is the ratio of Jg/J-total ’ 


it is now 


possible, knowing the surface temperature, to calculate the surface work 


function with greater accuracy than was previously possible. Because the 
body of the cathode was fabricated from quartz and the insert temperatures 
could be measured using an optical pyrometer, it was decided to conduct an 
experiment designed specifically to determine the work function of the 
electron emitting surface. If the emission mechanism is the one proposed, 
the surface work function which is calculated from the results should agree 


with accepted values for the materials used. 


The experiment consisted of first operating the cathode at a constant 


internal cathode pressure of '^5 torr and emission currents ranging from 0.8 
to 5.3 A, and then holding the emission current constant while operating 
the cathode at internal pressures ranging from 1.8 to 9.7 torr. The data 
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collected during the experiment were the current to each cathode surface, 
and the surface temperature of the insert segment which was doing most of 
the emitting. Using these data, the insert surface area, and the plasma 
properties obtained from previous tests, the surface work function (c|)g) of 
the insert segment was calculated from Eqs. 4, 6, 7, and 8. The results 
obtained from evaluation of these data are shown in Fig. 5 as the open 
symbols and dashed line, where the surface work function is plotted as a 
function of insert segment temperature. It should be noted that, although 
there were negligible temperature gradients in the axial direction of the 
rather short ('^2 mm long) insert segment, there were locally hot regions 
around the circumference which were as much as 56°C hotter than the average 
temperature used in calculating the work function. Because of these 
temperature gradients, another experiment was designed in which a smaller 
portion of the insert would be electrically and thermally isolated to 
provide data for calculating the work function. A sketch of such an insert 
with a test "patch" is shown in Fig. 6. The work functions calculated for 
this "patch" experiment are shown as the solid symbols and solid line in 
Fig. 5. 

The data of Fig. 5 suggest that the results for both the insert segment 
and the patch are similar and that in both cases the surface work function 
is a linearly increasing function of the surface temperature. Such a linear 
relationship between work function and surface temperature is, in fact, 
expected when the theoretical constant for A^(120 A/cm^°k^) is used to 
calculate (t)^ from the Schottky equation. In this case cj>^ is known as the 
"true" work function, as opposed to the "Richardson" work function, which 
is not a function of temperature but relies on an empirically determined 
value for the constant A . The relationship between these two work functions 
is discussed in Ref. 7. A representative value of the "true" work function. 



SURFACE WORK FUNCTION (eV) 
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Fig. 5. Effect of Temperature on Insert Work Function 





Fig. 6. 


Insert with Test Patch 
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(j)g, was not available for the tantalum insert coated with R-500 that was 
used in the experiment. However, Ref. 8 gives the work function of porous 
tungsten coated with a similar mixture of barium and strontium carbonates, 
as 


(f)^ = 1 .67 + 3.24 X 10"'+T(°C) . (9) 

This is plotted as the center-dash curve in Fig. 5. When the experimental 
results in Fig. 5 are compared to the center-dash curve their magnitudes 
appear to be reasonable (1.7 to 2.2 eV). Also both the experimental results 
and the center-dash curve display a linear increase with temperature though 
with different slopes. It might be pointed out that for the operating 
conditions in these tests the electric field adjacent to the emission 
surface, as estimated from Eq. 8, is ''^3X10^ V/cm. This means that when 
the cathode is operating, the field-enhancement results in an effective work 
function (jj^ which is reduced by '^^0.2 eV (Eq. 7) from the surface work 
function (p^ shown in Fig. 5. 

It should be noted that an incorrect estimate of the emission surface 
area can cause a change in the slope of the linear work function-temperature 
plot like the one shown in Figure 5. If the extent of the emission surface 
area is overestimated then the slope of the work function curve will be 
increased over what it would be if the true emission area were used for 
calculating the work function. The difference in slope is given by the 
expression 


m 


act 


„ _ k . /_e^\ 

"'est " e 

act 


( 10 ) 


where m^^^ and m^^^ are the actual and estimated slopes of the curves, and 
^act ^est corresponding emission areas. To account for the 

difference in slope between the solid curve in Fig. 5 for the patch data 
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and the curve based on Ref. 8 solely on the basis of having overestimated 
the emission area would imply that the estimated area was 2300 times the 
actual emission area. Since the patch temperature was quite uniform it is 
rather unlikely that the emission was coming from such a small area unless 
the work function was highly nonuniform across the surface. A more likely 
explanation is that the actual surface work function for tantalum coated 
with R-500 deviates somewhat from that of porous tungsten coated with 
barium and strontium carbonates. 

Future Work 

The results shown in Fig. 5 indicate that the emission process is 
characterized by a work function that is reasonable because of its magnitude 
and its linear dependence on temperature. However, it is apparent from the 
discussion above that the exponential sensitivity of the emission current 
to the work function makes this evidence insufficient to permit one to 
determine conclusively whether field-enhanced, thermionic emission is the 
only emission mechanism acting. Such a judgement could be made, however, 
if the actual work function were accurately known for the emitting surface. 
One method of determining this work function would be to make an inde- 
pendent measurement of the work function of the surface without the discharge 
present. Because of the extreme sensitivity of the surface work function 
to contamination and operating history, it is considered imperative to 
make this measurement without having to remove the insert from the cathode. 

An experiment to perform such an in-situ measurement has been proposed and 
is currently in progress. The basic features of the apparatus to be used 
in this test are illustrated schematically in Fig. 7. The cathode would 
first be operated in a normal manner and the emission current and temperature 
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Fig. 7. Apparatus for In-situ Measurement of Insert Work Function 
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of the insert "patch" would be measured. An ion collecting probe would 
be used to measure plasma density and electron temperature adjacent to 
the "patch". The results of this test would be used as input to the 
field-enhanced, thermionic emission model, and the work function associated 
with this operating mode would be calculated. Next the discharge and 
mercury flow would be turned off, and the external heater would be used 
to heat the insert and "patch" to emission temperatures. The collector 
anode shown in Fig. 7 would then be biased positively to collect the 
emission current. The measured patch current and temperature would now be 
used to calculate a second, independently determined work function for the 
"Patch". Comparison of the two work functions should provide a conclusive 
test as to whether or not field-enhanced, thermionic emission is the dominant 
emission process in the hollow cathode. 

As mentioned earlier, the model lacks the ability to explain the 
location and extent of the emission region or to make a priori predictions 
of the plasma properties. Additional efforts are planned to model the ion 
production region in greater detail in order to fill in these gaps. This 
work will look at modeling the ion production and loss rates within the 
ion production region using the present model to specify the boundary 
conditions for the region. It is hoped that this will allow prediction of 
the plasma properties in this region and provide some insight into a 
criteria for determining the emission location and extent. 

Concl usions 

The proposed phenomenological model describing orificed, hollow cathode 
operation is in good agreement both qualitatively and quantitatively with 
a number of recent experiments. The model's description of the ion pro- 
duction region for the internal discharge as the region bounded by the 
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electron emission surface is supported by the experimental results. 

Further, the model makes reasonable and self-consistent predictions of 
the plasma density in that region. The current to the cathode emission 
surface makes up 85% of the total emission current with about 20% 
of that being attributable to the ion current to the emission surface. 

The remaining 15% of the total emission current is due to ion currents to 
various cathode and facility ground surfaces. These surfaces do not 
collect a significant electron current, nor do they appear to produce any 
appreciable electron emission. 

The "true" surface work function of the emission surface calculated 
from experimental data using the model shows the expected linear dependence 
with temperature and is of a reasonable magnitude. Though this is not 
sufficient to prove conclusively that field-enhanced, thermionic emission 
is the dominant emission mechanism in this type of cathode, it does tend 
to support this assumption. An additional experiment which could possibly 
provide such proof is suggested. 

The orifice region sustains potential drops of a few volts and appears 
to be an ion production region. The plasma density in this region shows a 
different sensitivity to internal pressure changes than does the internal 
discharge, making it likely that this region is isolated from the internal 
discharge by phenomena such as a double-sheath or gas dynamic effects. 

Finally, though the proposed model falls short of making a priori 
predictions of plasma properties within the cathode, it can predict the 
emission surface temperature if the plasma potential, electron temperature, 
emission surface area and surface work functjon are assumed or specified. 
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HOLLOW CATHODE PLASMA OSCILLATIONS 
Dan Siefried 

Work directed at refining the Langmuir probe technique used for 
making plasma measurements inside the hollow cathode identified plasma 
oscillations in the megahertz range that appear to be an inherent feature 
of orificed, hollow cathode operation. Initially investigated to de- 
termine their effect on Langmuir probe results, these oscillations were 
found to be an interesting and possibly important phenomena in their own 
right. Although these oscillations were not studied in detail, it was 
determined that they are due to some plasma phenomenon occurring in the 
orifice or downstream of the orifice and that they are not dependent on 
power supply characteristics. 

Experimental Results 

The quartz test cathode used in the experiments to be discussed below 
was similar to the one used in earlier experiments and is shown in Fig. 8. 

The cathode was operated in the normal manner while various parameters such 
as discharge current, mass flow rate and orifice diameter were varied and 
parameters such as anode voltage, keeper voltage and probe floating potential 
were displayed on an oscilloscope. The frequency and amplitude of the 
observed oscillations were found to be affected to some degree by each of 
the parameters that was varied and also by operating mode (spot or plume). 
The oscillations varied from very well defined, coherent waveforms at 
frequencies around one megahertz to what appeared to be either random noise, 
or possibly oscillations at frequencies beyond the range of the oscilloscope 
which was being used ('v50 megahertz). These oscillations were not investi- 
gated in sufficient detail to define the effects of various parameters on 
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the oscillations, although it was noted that the most well-defined wave- 
forms were found when the cathode was operating in the spot mode at currents 
just above the transition point between the plume and spot modes. A 
photograph of the oscilloscope trace displaying a typical waveform under 
these conditions is shown in Fig. 9. This photo was taken for a test 
cathode having a 0.76 mm diameter orifice operating at a total emission 
current (I^) of 3.5A, a discharge voltage (V^^) of 16.6v, and a mass flow 
rate (m) of '\'100 mA. The oscillations have a frequency of '\^0.9 MHz and an 
amplitude of '\^0.15v. The oscil lations in this case were sensed with a 0.8 mm 
diameter spherical probe with the probe at floating potential ('v^2.5v) and located 
'v^lO mm upstream of the orifice plate. The oscillations were also present 
on the keeper and the anode at the same frequency but with larger amplitudes. 
For comparison, Fig. 10 shows the oscillations present in the plume mode 
at the same flow rate for a total emission current of 2.3 A and a discharge 
voltage of 21.8 v. In this case the waveform is not as well defined and 
appears to have a frequency of 0.67 MHz and an amplitude of 'vO.06 v (note 
that the vertical or voltage scale in this figure is 0.1 v/division as 
compared to 0.05 v/division in Fig. 9). The probe voltage fluctuations 
at a higher current (l^ = 5.3 A; = 17.5 v; spot mode) are shown in 
Fig. 11. Results shown here were typical of higher current conditions 
where the coherent oscillations appeared to break down into random noise. 

Note that the horizontal, or time scale in this oscillograph is different 
than that in the previous ones. However, even when displayed on the same 
time scale, that is 0.5 ps/division (the limit for this oscilloscope), no 
coherent oscillations could be detected at this higher current condition. 

Since the test cathode was made of quartz with a tantalum orifice 
plate, it was possible to disconnect the keeper and normal anode and 
collect the discharge current at the upstream side of the orifice plate 
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Fig. 9. Probe Floating Potential _ FI actuations 

(I^ = 3.5A, = 16.6v, m = 100mA, Spot Mode) 
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Fig. 10. Probe Floating Potential . FI actuations 

(I^ = 2.3A, = 21 .8v, m = 100mA, Plume Mode) 
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Fig. 11. Probe Floating Potential Fluctuations 

(I^ = 5.3A, = 17. 5v, m = 100mA, Spot Mode) 



Fig. 12. Probe Floating Potential -El ectron Collection at Orifice Plate 
(I^ = 2.3A, = 12.lv, m = 100mA, Plume Mode) 
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by utilizing it as an anode. When this was done it was found that there 
was no visible discharge present downstream of the orifice. In this con- 
figuration the plasma oscillations disappeared completely and the internal 
discharge was quiescent. This is illustrated graphically by the oscillo- 
graph shown in Fig. 12 which shows no A.C. component for the probe 
floating potential when the oscilloscope was on its most sensitive scale. 

When the keeper was reconnected and its voltage increased, a portion of 
the discharge again coupled through the orifice to the keeper and the 

oscillations reappeared, 
g 

Byers, studying the effect of power supply impedance on the SERT II 
neutralizer, also found coherent oscillations present at similar frequencies. 
Byers found that the oscillations were influenced by the output impedance 
of the cathode discharge power supply, depending in particular whether this 
impedance was capacitive or inductive. To determine if the oscillations 
observed in the present experiment were either caused or affected by power 
supply characteristics a test was run in which battery packs were substituted 
for the regulated power supplies normally used. In this test the cathode 
was started and operated in the normal manner and the oscillations were 
displayed on an oscilloscope. Both the keeper supply and the anode supply 
were then replaced with battery packs set to provide the same steady state 
operating conditions as those realized with the conventional power supplies. 
All peripheral supplies (heaters, vaporizer, etc.) were then shut off 
momentarily. The plasma oscillations were not affected either in amplitude, 
frequency, or waveform by this use of battery power supplies. 

Discussion 

These tests indicate that the oscillations are not due to power supply 
interactions and are apparently being produced either within the orifice 
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region or downstream of the orifice in the keeper region. One possible 
source of oscillations in this frequency range is an ion acoustic wave. 
Crawford^’ in an investigation of the double sheath which forms at a 
discharge constriction has identified oscillations related to standing, 
radial , ion acoustic waves. The expression for the fundamental frequency 
for this type of phenomenon in a cylindrical geometry is 


^0 ^ ird m, ‘ 


01 ) 


where d is the diameter of the constriction, y is the adiabatic com- 
pression coefficient (either 1, 5/3, or 3 but probably 3 for this case), 
k is the Boltzmann constant, is the electron temperature, and m^. is the 
ionic mass. For mercury under conditions where y = 3 and the electron 
temperature is 0.8 eV 


^0 = 


820 


( 12 ) 


If such a standing wave is present in the orifice, then for an orifice 
diameter of 0.76 mm the fundamental frequency for the oscillation would 
be f^ = 1.08 MHz which is in good agreement with the observed frequency. 

This still leaves a question regarding the mechanism for exciting 
these oscillations. Crawford suggests that the double sheath which he 
observed at a plasma discharge constriction may be set into motion, like 
a vibrating membrane, by the electrons which are being accelerated through 
it. Experiments discussed in the preceding section showed that the 
orifice region is probably an ion production region and suggested that it 
may indeed be separated from the internal discharge by such a double sheath. 
One way of viewing this excitation mechanism is that as the double sheath 
is formed electrons are accelerated across the sheath with sufficient 
energy to ionize neutrals in the orifice region. This establishes a higher 
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plasma density downstream of the sheath in the orifice constriction - 
such potential and plasma density differences between two regions are 
characteristic of the formation of a double-sheath boundary between two 
plasmas. The oscillations are then excited as the ions in the orifice 
region diffuse to the orifice wall and upstream through the sheath, thereby 
causing the double sheath to collapse or at least for the potential across 
the sheath to fluctuate. If this is the case, the characteristic time for 
an ion to escape from the production region should be related to the 
frequency of the oscillations. Based on the Bohm velocity for an electron 
temperature of 0.8 eV and using the orifice diameter as. a characteristic 
length, a mean transit time for an ion formed in this region to be lost to 
the wall or upstream through the sheath can be calculated and is found to 
be on the order of 1.22 ysec. This would correspond to a frequency of 
0.82 MHz which again is in reasonable agreement with the frequency of the 
observed oscillations. This calculation is presented for its physical 
descriptiveness and is just a simplification of Eq. 11. Equation 11 is in 
fact the ion acoustic velocity divided by the characteristic dimension 
with constants dependent on the geometry and boundary conditions. 

Concl usions 

Plasma oscillations appear to be an inherent feature of orificed, 
hollow cathode operation. The oscillations are due to some plasma process 
taking place either in the orifice region or downstream of this region. 

The oscillations are not due to power supply interaction with the cathode. 
The oscillations may be in the form of standing, radial , ion acoustic 
waves excited by fluctuations in the potential across a double sheath as 
ions are produced and then lost to the walls in the orifice region. 
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PROBE MEASUREMENTS INSIDE THE ORIFICE[; HOLLOW CATHODE 

Dan Siegried 

Considerable effort was directed at developing and verifying a probe 
technique and analysis suitable for determining plasma properties within 

3 

the orificed, hollow cathode. Previously, measurements have been made 
within the cathode using a cylindrical Langmuir probe with a very fine 
0.076 mm diameter tungsten electrode. This gave reasonably good results 
in regions where the plasma density was 10^‘*cm“3 or less. However, in 
regions where the emission is taking place which are of greatest interest 
the plasma densities are greater than lO^^^cm"^ and the very fine wire probe 
is found to burn up when biased near plasma potential. In order to alleviate 
this problem, a new technique using a spherical probe having a larger sur- 
face area and an analysis based on the ion saturation region of the probe 
trace was investigated. 

Because barium oxide compounds are used to reduce the surface work 
function in the cathode, probe contamination was also an important concern. 
Probe circuitry was developed which incorporated a probe cleaning bias and 
allowed a rapid sweep of the probe potential. This was used as a tool to 
examine the effects of contamination on the probe results. The effect of 
plasma oscillations on the probe results was also investigated. 

Spherical Probe Design 

A number of probe designs were tried with varying degrees of success 
before settling on the one shown in Fig. 13. In this design, the probe 
electrode consists of a ^0.8 mm diameter sphere formed by melting the end 
of a 0.25 mm diameter tungsten wire in a D.C. helium arc. The 2 mm 
diameter quartz insulator is drawn down to a diameter of '\'0.4 mm over a 
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Fig. 13. Spherical Probe Design 
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length of about 2 cm to minimize its perturbing effect on the plasma. 

The end of the insulator is wrapped with a couple of layers of 0.025 mm 
tantalum foil. This acts as both a thermal shield and a sputter shield 
for the end of the insulator. The wire is kinked just behind the sphere 
to give a tight fit in the insulator. This helps to minimize the change 
in probe surface area as the probe lead heats up and elongates during 
probe operation. When used for making measurements inside of the hollow 
cathode, the probe is inserted through a feedthrough on the upstream end 
of the cathode mounting fixture in such a way that it can be moved to 
various positions along the axis of the cathode. Figure 13 is approximately 
to scale and shows the insert diameter in correct proportion to the probe 
dimensions . 

Data Analysis 

Earlier probe measurements made inside the hollow cathode required 
biasing the probe to potentials above the plasma potential to obtain a 
probe trace such as the one shown in Fig. 14. As mentioned in the intro- 
duction, biasing the probe to plasma potential in regions where the plasma 
density is very high causes the probe to burn up. The large surface area 
of the spherical probe (Fig. 13) however makes it sensitive to the ion 
saturation current over the wide range of plasma densities found inside 
the cathode. It also provides sufficient information to facilitate 
determination of the plasma properties by biasing the probe to the floating 
potential (see detail of Fig. 14). Operation at or below floating potential 
has the additional advantage of causing less of a perturbation to the plasma 
than does the large electron current drawn to a probe at plasma potential. 
The use of the ion saturation region of the probe current-voltage character- 
istic requires a different method of analysis than does the usual Langmuir 
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probe trace. A method for performing this analysis will be presented, but 
first it will be useful to discuss some considerations which affect the 
type of analysis used. 

In the past, characteristics taken from probes inside of the hollow 
cathode have been analyzed based on Langmuir's theory. This analysis 
applies in the limit where the sheath is collisionless and is thin with 
respect to probe size. This implies a low pressure, slightly ionized plasma. 
The hollow cathode plasma, however, is a relatively high pressure (a few 
torr) plasma with a high plasma density (10^^ to 10^^ cm"^). Even though 
the pressure is high, it has been assumed in the past that because of the 
very thin Debye length (due to the high plasma density), the normal analysis 
could be applied with reasonable results. It still appears that this is 
the case, based on the argument contained in the following paragraphs. 

The conditions found in the hollow cathode are such that, although 

the sheath is thin compared to the probe dimensions, the high neutral 

density within the cathode results in particle mean free paths which are of 

the same order as the probe dimensions . This usually causes a presheath to be 

established between the bulk plasma and the sheath itself in which particles 

approaching the probe undergo numerous collisions. Significant gradients 

in the plasma density and potential would normally be expected in the 

presheath as a consequence of these collisions. When the presheath is 

present, the conventional Langmuir analysis is usually inaccurate in that 

it significantly underestimates both the plasma potential and the plasma 

density of the undisturbed bulk plasma. A survey of literature has shown 

11 12 

that there are a number of theories ’ which are available for analyzing 
the probe data in a situation where the sheath is thin and the high neutral 
density causes a significant presheath. In the case of the hollow cathode. 
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however, the region surrounding the probe is also one in which the ion 
production rate is large compared to the rate of ion loss to the probe. 
Further, the probing technique being proposed here involves the use of 
only the ion saturation portion of the characteristic probe curve so that 
large electron currents are not drawn to the probe. Under these conditions, 
it is believed that the presheath plasma gradients discussed above cannot 
be established. 

Under typical hollow cathode operating conditions it is estimated that 
the production rate of ions in a volume of the order of the presheath volume 
is considerably greater than the loss rate of ions to the probe. Where this 
is the case it is considered justifiable to assume that no significant plasma 
gradients can be sustained and the conventional thin sheath analysis can be 
applied. This is almost certainly the situation in the high plasma density 
{10^3 to 10^^ cm"^) region of the cathode where electron emission and ion 
production are concentrated. Upstream of this region where more moderate 
(10^^ to 10^^ cm"^) plasma densities are expected some intermediate situation 
probably occurs in which the ion production rate is not great enough to 
compensate for the plasma gradients set up as ions diffuse to the probe. 
However, in order to standardize and simplify the analysis,' and because 
the method of approach in the intermediate situation is uncertain, it is 
proposed that all of the probe results obtained inside of the hollow cathode 
be analyzed using the relatively simple, thin sheath, collisionless type of 
analysis . 

Given the ion saturation portion of the probe characteristic shown in 
the detail of Fig. 14, the analysis for the thin sheath probe in an un- 
preturbed plasma (i.e. no presheath) is quite simple. First, the electron 
temperature is determined in the standard manner from the exponential portion 
of the trace. This is done by making a semilog plot of the electron current. 
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ig (shown in the detail of Fig. 14) above the extrapolated ion saturation 
current, against probe potential. The electron temperature is then pro- 
portional to the inverse slope of this line. The plasma density is 
determined from the ion saturation current (i^.) based on the assumption 
that the ions enter the sheath with the Bohm velocity. This can be ex- 
pressed as 
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where Ap is the probe surface area, e is the electronic charge, k is 
Boltzmann's constant, T^ the electron temperature, and m^. is the ion mass. 


The parameter 


is dependent on the ratio of ion temperature to electron 


temperature and the size of the probe in relation to the extent of the plasma 

1 3 

but usually has a value of '\^0.6. The plasma potential can be estimated 

from the floating potential based on the condition that the ion and electron 

currents to the probe must be equal when the probe is floating. This leads 

1 3 

to the following expression for the plasma potential: 
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where is the floating potential, m^ the electron mass, and the rest of 


the symbols are as previously defined. 


Probe Contamination 

A common problem with probe traces taken in an oxide cathode environment 

1 4 

is that the probe will be contaminated by the low work function compounds. 
This can cause two problems: 1) the formation of resistive oxide coatings, 

and 2) a reduction in the probe work function. The resistive coatings 
cause a change in the shape of the characteristic, and thereby, give a false 



47 


indication of the velocity distribution of the electrons in the plasma. 

The probe work function affects the position of the probe characteristic 
with respect to the reference potential , and can therefore affect the 
value determined for the plasma potential. This comes about because a 
probe characteristic must be recorded with reference to some other elec- 
trode exposed to the plasma, and the work function acts essentially as a 
contact potential between the electrode and the plasma. Ideally the 
reference electrode would have the same work function as the probe. 
However, if this is not the case then the probe characteristic will under- 
go a parallel shift along the voltage axis depending on the difference in 
work functions between the probe and the reference electrode. Generally, 
this parallel shift is small compared to the plasma potential and is of 
1 ittl e consequence. 

A more important problem occurs if the probe work function changes 
during the period while the trace is being taken. This will result in a 
distortion of the shape of the probe characteristic again leading to false 
indications of the electron velocity distribution. These problems are 
generally addressed by taking precautions to maintain the probe in a clean 
condition and then recording the probe trace rapidly enough so that the 
probe surface conditions do not vary significantly. In the present case, 
this was accomplished by designing and building probe circuitry that 
maintains the probe at a bias potential sufficient to clean the probe and 
then allows one to sweep the probe potential rapidly to obtain the probe 
characteristic. The circuitry used to do this is described in Appendix A. 
This probe circuitry has proven to be a valuable tool for producing clean, 
consistent probe traces and for investigating the effects of probe contam- 


ination. 
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Using the probe circuitry described above, it was determined that 
probe contamination probably had been affecting probe traces taken inside 
of the cathode. Numerous probe traces were recorded with sweep times of 
''^1 sec. (near the response limit of the recorder) while various cleaning 
biases were applied to the probe. The most effective cleaning was ob- 
tained by maintaining the probe about a hundred volts negative of cathode 
potential when a trace was not being recorded. This assured ion energies 
sufficient to sputter clean the probe surface. This was found to be 
preferable to electron heating of the probe to incandescence which was 
found to alter plasma conditions significantly and is generally not 
effective in removing resistive oxide coatings because of their very high 
evaporation temperatures. The effect of contamination on the probe was 
detected by varying the sweep time and also by varying the time delay 
between switching off the cleaning bias and starting the trace sweep. If 
the probe trace was swept in a few seconds immediately after turning off 
the probe cleaning circuit the traces were quite repeatable. However, if 
the trace sweep was taken slowly (5-10 seconds) or if the sweep was delayed 
for times on the order of 5-10 seconds, the floating potential of the trace 
would shift and the shape of the trace would begin to change. This 
suggested that in the cathode environment the time constant for probe con- 
tamination is on the order of 10 seconds and that the probe will probably 
not be significantly contaminated if the trace is swept within a few 
seconds immediately after the cleaning circuit is turned off. 

At this point it was concluded that by using the above technique the 
probe should be substantially free of contamination. However, the probe 
traces that were being recorded indicated electron energy distributions 
that appeared to be non-Maxwel 1 ian at least in the lower energies (i.e.. 
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near plasma potential). As mentioned earlier, one possible cause for 
this is a change in probe work function as the trace is swept. If some 
type of contamination affecting the work function were still present, 
then a distorted probe trace would be obtained as the probe was heated 
causing the contaminants to evaporate during the trace sweep. To check 
whether or not this was occurring, it was considered desirable to run an 
experiment in which probe traces obtained at increasingly fast sweep speeds 
were compared. It was decided that the sweep speed should be increased up 
to the point where no significant probe heating could be taking place. 

Probe traces obtained more rapidly than second could not be obtained 
using the x-y recorder, so the traces were taken using a storage oscillo- 
scope. It was found that no difference in the probe traces could be de- 
tected for sweep times ranging from a few seconds down to 2 milliseconds 
(the response limit of the differential amplifier used for measuring the 
probe current). Probe temperature changes were estimated based on calcula- 
tions of the thermal response of the probe. Under the conditions of the 
experiment, it was estimated that the probe temperature increase was a 
maximum of '^^20°C at 1 second and only a few tenths of a degree for a 
millisecond sweep time. This indicates that under the proposed probing 
technique the probe is probably in a clean condition and remains in that 
condition without significant work function changes during the few seconds 
required for recording the trace. The indicated non-Maxwel 1 ian nature of 
the electron energy distribution is apparently not attributable to con- 
tamination effects. 

Plasma Oscillations 

As discussed in the preceding section of this report, plasma 
oscillations in the megahertz range are an inherent feature of orificed. 
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hollow cathode operation. A question was raised as to whether or not 
these oscillations could be causing the indication of a non-Maxwell ian 
electron distribution. In order to examine this question, an analysis 
was undertaken which looked at a Taylor series expansion of the electron 
current collected by a probe under the influence of a sinusoidally varying 
voltage. This analysis showed that, if the electron population is truly 
Maxwellian, the probe traces will have the same characteristic shape and 
indicate the same electron temperature regardless of whether or not 
oscillations are present. However, the oscillations will cause the elec- 
tron density indicated by the probe trace to be greater than the actual 
density by a factor of 
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where is the amplitude of the plasma oscillation. Under prevailing 


conditions in the cathode, however, this effect results in a density 
variation of only a few percent. The analysis suggests that either the 
plasma electron population is truly non-Maxwel 1 ian or that there is another 
cause for the non-Maxwel 1 ian probe traces. The actual oscillations 
observed in cathodes have been found to be more saw tooth shaped than 
sinusoidal, but the basic results of this analysis are probably not 
changed significantly by this fact. 


Additional Considerations 

A few additional comments are in order regarding the non-Maxwell ian 
nature of the probe traces. For most regions of the cathode the collection 
of the energetic electrons, that is the portion of the electron population 
collected at probe potentials between cathode potential and floating 
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potential, does in fact follow to a reasonable extent an exponential or 

Maxwellian pattern. It is only at lower electron energies, that is for 

probe potentials between floating potential and plasma potential that the 

11 12 

electron distribution appears to be non-Maxwell ian. Collisional effects ’ 

are known to distort the exponential nature of the probe traces at these 

lower energies, and this was in fact the impetus for examining the 

collisional effects that were discussed in the section on probe analysis. 

However, even when collisional effects are present the electron retarding 

11 12 

region near ion saturation is expected to show a Maxwellian behavior ’ 
and it is this region which is used for determining the electron temperature. 
This suggests that the proposed analysis using the exponential portion of 
the probe trace found at potentials below the floating potential where 
probe perturbation effects are smallest is probably justified and gives the 
best estimate for the electron temperature. However, this statement should 
be qualified by acknowledging what it implies physically. The electrons 
being collected at probe potentials near cathode potential make up only a 
small fraction of the total electron population. Moreover, if the cathode 
model proposed in the preceding section is correct, then the electrons 
produced by acceleration from the emission surface into the plasma would be 
expected to form a reasonably mono-energetic population with an energy in 
the range that would be detected by a probe near cathode potential. This 
mono-energetic population is estimated to be as much as a few percent of 
the total electron population, so that based on the present physical under- 
standing of the cathode, one would predict that the probe trace, rather than 
being Maxwellian near cathode potential , would be non-Maxwel 1 ian probably 
showing a Maxwellian plus mono-energetic type of behavior. However, there 
has been no definite evidence that this is the case. At this point, it 
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appears that the best that can be done is to determine the electron 
temperature based on the energetic electrons collected at potentials 
below floating potential, and to assume that the bulk of the electron 
population is in fact Maxwellian at that temperature. 

Concl usions 

A spherical probe design suitable for making measurements inside the 
hollow cathode has been presented. The proposed method of analyzing the 
probe characteristics is shown to be reasonable for the plasma conditions 
prevailing within the cathode. Probe circuitry and operating procedure 
are presented which are capable of maintaining the probe in a clean 
condition between traces and of minimizing the effect of contamination 
during the recording of a trace. Plasma oscillations present during 
cathode operation do not have any substantial effect on the electron 
temperature or plasma density indicated by the probe characteristic. 
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SPUTTERING RATES OF DISCHARGE CHAMBER MATERIALS 
Susumu Masaki 

The lifetimes of discharge chamber components are frequently de- 
termined by the erosion rates of those components. This erosion is 
caused by discharge chamber ions (primarily doubly-charged ones) which 
strike these surfaces and sputter material away. It has been observed 
that impurity atoms present in the discharge chamber and on discharge 
chamber surfaces can affect the erosion rates of these components. Nitrogen 
for example has been shown to reduce sputtering rates of discharge chamber 

1 5 

materials greatly even when it is present at very small partial pressures. 

It is considered likely that even the mercury propellant may alter the 
erosion rate of a surface from that which might be expected for the case of 
doubly-charged mercury ions bombarding a clean surface. The measurement 
of sputter erosion rates has become a matter of considerable interest for 
these reasons. 

Sputtering rates can be determined in a discharge chamber environment 

by actually operating a discharge chamber for a known period of time, 

measuring the depth of sputter erosion and then calculating the erosion 

rate. In order to obtain accurate values however very long test times are 

required. In order to facilitate shorter tests, other less direct methods 

have been pursued. One of these methods involves installing multilayer 

1 6 

erosion test patches. These patches are made by alternating very thin 
layers (several hundred angstroms thick) of the material being tested with 
copper layers of similar thickness. By counting the number of layers of 
known thickness that are eroded during a few hours of testing the erosion 
rate can be determined. A question can however be raised about the validity 
of such a test, namely; can the erosion rates of copper and the test 
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material in the layered configuration be assumed to be the same as those 
determined from measurements on pure copper and pure test metal? If this 
can be assumed then the test can be used, but if the presence of the copper 
induces an effect such as coning^ ^ then the results may not be accurate. 

In order to determine the sputtering yields of various thruster 
materials at the ion energies observed in discharge chamber plasmas, and 
to investigate the effects of such impurities as copper and nitrogen on 
these erosion rates the tests described in the following paragraphs were 
conducted . 

Apparatus and Procedure 

These tests were conducted in an 18 in. dia. bell jar facility. The 

ion source used was a 15 cm dia. multipole design with its grid masked down 

to a 10 cm dia. This source produced current density profiles like those 

shown in Figure 15. Test samples to be sputter eroded in the beam were 

placed 6 cm downstream of the grids within the central 3 cm radius circle 

where the current density was very uniform. The source used either argon 

or mercury working fluids and the argon or mercury pressures were typically 
-4 5 

'\^5X10 or 3X10 torr respectively during operation. Figure 16 shows 
the ion source with the test specimens. The specimens consisted of a flat 
glass substrate onto which a uniform layer of test material several thousand 
angstroms thick had been sputter deposited. These specimens were fabricated 
locally and considerable care was exercised to insure that they were kept 
clean. In order to minimize oxidation they were kept under vacuum except 
when they were being transferred between vacuum facilities. After fabrica- 
tion the sputter coated substrate was partically covered by two 1 mm 
thick graphite shields that had been undercut in the manner suggested in 
Figure 16. As a test specimen was exposed to the ion beam, material was 
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Fig. 15. Current Density Profile 
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Fig. 16. Sputter Erosion Test Apparatus 
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protected where it was covered by the graphite shield and eroded where it 
was not. As the figure suggests ail test samples were exposed to normal 
incidence ion beams. This resulted in a gradual step having a height which 
could be measured using an optical interferometer (Angstrometer) after 
exposure to the beam had been terminated. Test specimens were 1 cm by 
4 cm and five specimens could be fitted into the uniform current density 
portion of the beam so they could be exposed simultaneously. Several kinds 
of glass substrates ranging from optical flats to photographic plate glass 
were tried during the course of the tests and photographic plate glass was 
found to be sufficiently flat so non-uniformities in the glass did not cause 
errors in thickness measurements made with the optical interferometer. 

In conducting the test, the test samples were installed in the vacuum 
facility and the system was pumped down to about 1 X 10" torr. The source 
was started with the shutter shown in Figure 16 covering the samples, and 
stable source operation at the desired current density and ion energy was 
established. The shutter was then opened and the test was initiated. 
Typically the samples were exposed to an ion beam of constant current density 
and energy for about 100 minutes. During exposure the samples were cooled 
radiatively and their temperature rose rapidly to a value of 120 to 140°C 
when they were exposed to the ion beam. Typically their temperature had 
risen through 90% of the final temperature rise within about 7 minutes 
after the shutter had been opened. The samples were allowed to cool for 
several hours after they had been sputtered. They were then removed from 
the vacuum facility and the height of the steps machined into them at the 
edges of the graphite shields were measured. 

To establish one data point, five separate measurements of the step 
height were made on each of two separate machined steps. The mean value of 
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the step height was then computed and this was divided by the test duration 
to determine the mean erosion rate. Typically the standard deviation 
associated with these measurements was i 10% or less. Using this mean 
erosion rate and the beam ion current density at the location of the test 
sample, the sputtering yields associated with each test specimen were 
computed. Because of the low net acceleration voltage and large resultant 
beam divergence associated with the tests, considerable care had to be 
exercised to assure that current densities used as input to this computation 
were correct. This was accomplished by correlating shutter-probe current 
density measurements with those obtained at the sample location at each 
ion energy condition. Corrections to the current density to account for 
fast moving charge-exchange atoms which should also cause sputtering erosion 
were also incorporated. 

The singly-charged ion energies for these tests were determined by the 
ion source anode-to-ground electrical potential difference. For these tests, 
values of this potential difference that were twice the anode-to-cathode 
potential difference observed in thrusters were selected. Under this con- 
dition singly-charged ions in the beam of the test set up have an energy 
equal to the doubly-charged ions that are the dominant cause of sputter 
erosion in thruster discharge chambers. Doubly-charged ions in the test 
set up of Figure 16 would have twice this energy and would introduce an 
error in the sputtering yield results. The discharge voltage was held 
at the lowest value where a stable discharge would still occur (35v for 
mercury and 45v for argon) to minimize this error. 

In a second series of tests nitrogen was bled into the vacuum facility 
to increase the pressure by a specified amount after thruster operation had 
been stabilized. The sputtering of each test sample in this environment was 
measured using the same procedure described in the preceding paragraphs. 
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By varying the change in vacuum chamber pressure caused by the nitrogen 
(nitrogen partial pressure) the effect of this parameter on the sputtering 
yields of various thruster materials was investigated. The tests were 
conducted for one ion energy condition when operating with each working 
fluid (105 eV for mercury and 115 eV for argon). 

In a final series of tests, samples were fabricated using alternated 

O 

500 A layers of copper and molybdenum and of copper and 304 stainless 
steel. These test patches were exposed to argon ion beams ranging in 
energy from 75 to 115 eV at various nitrogen partial pressures. Molybdenum/ 
copper layered test samples were also exposed to 105 eV mercury ion beams 
at various nitrogen partial pressures. After the samples had been sputtered 
the erosion depth was measured using the aforementioned procedure. 

Results 

Test results showing the effect of argon ion energy in the range 75 

to 135 eV on the sputtering yields of tantalum, molybdenum, 304 stainless 

steel and copper are given in Figure 17. Figure 18 shows the effect of 

mercury ion energy on the sputtering yields of the same four materials. 

The argon ion sputtering yield data of Figure 17 for copper are in reasonably 

1 8 

good agreement with those obtained by Wehner, the yields for tantalum are 
about half of the values measured by Wehner^^ and Stuart and Wehner^ ^ at 
135 eV and are in agreement near 75 eV. For mercury the results of 

1 g 

Figure 18 can be compared with sputtering yield data obtained by Wehner 

20 

and Askerov and Sena. The results of these researchers are in reasonably 
good agreement with each other while the sputtering yields of Figure 18 are 
in general higher than theirs. The reasons for these differences are not 
apparent, but one possible explanation for the higher erosion rates observed 
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Fig. 18. Mercury Ion Sputtering Yields 
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in the present mercury experiment is that doubly-charged ion erosion was 

not negligible as has been assumed. It is interesting to note that the 

second derivative of sputtering yield with mercury ion energy seems to be 

positive while for argon it seems to be negative in both the present data 

1 8 

and the results of Wehner. 

The effect of nitrogen partial pressure on the sputter yields of the 

various thruster materials being investigated at argon and mercury ion 

energies of 115 and 105 eV respectively is shown in Figures 19 and 20. 

Each data point shown is determined by averaging four separate measurements 

of the erosion depths on each of two separate steps machined by the ion 

beam. The scatter in these data about the mean values given in Figures 19 

and 20 suggest an error of t 20% in the data is likely. The data show the 

expected tendency for the sputtering yields to drop off substantially as 

the partial pressure of nitrogen is increased. 

The data of Figures 19 and 20 are open to question at the low nitrogen 

partial pressure levels because the ultimate pressure in the facility, 

—6 

where the sputtering was effected, was 'v 1 X lO” torr. Other investi- 
15 21 

gators ’ have shown using comparative spectroscopic measurements that 
the onset of the effects of air and nitrogen on sputtering yields occur 
at pressure at or below 1 X 10" torr. The data given here suggest that an 
approximate halving of the sputtering yeilds can be expected for most of 
these materials as the nitrogen partial pressure is increased from <1 X 10" 
to a few tens of millitorr. 

Erosion rate data obtained from test samples fabricated using alternate 
layers of copper and molybdenum or 304 stainless steel and copper are shown 
in Table V. Columns two and three of this table show the nitrogen partial 
pressures and argon ion energies at which sputtering was effected. The 
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TABLE V Sputter Erosion Characteristics 
of Multilayered Test Samples in an Argon Ion Beam 


Test 

Nitrogen 

Argon Ion 

Erosion 

Test 

Calculated 

Sampl e 

Pressure 

Energy 

Deot 

Duration 

Erosion Time 

Configuration 

(torr) 

(eVj 

0 

(A) 

(min) 

(min) 

500 A 



5930* 



304 SS/Copper 
Layers 

<6 X 10“® 

115 

(3020/2910) 

100 

96 

500 A 

A A 


2945 



Moly/Copper 

Layers 

<6 X 10’^ 

75 

(1500/1445) 

120 

129 

500 A 



4128 



Mol y/ Copper 
Layers 

<6 X 10"^ 

95 

(2128/2000) 

1100 

92 

500 A 

<6 X 10'® 


4551 



Moly/Copper 

Layers 

115 

(2500/2051 ) 

100 

105 


500 A 



4172 



Moly/Copper 

Layers 

5 X 10"'^ 

115 

(2172/2000) 

100 

106 

500 A 



2598 



Moly/Copper 

Layers 

1 X 10"^ 

115 

(1 500/1098) 

100 

113 

500 A 



1841 



Moly/Copper 

Layers 

2X10"^ 

115 

(1000/841 ) 

1000 

92 

500 A 



1473 



Moly/Copper 

Layers 

3X10^ 

115 

(973/500) 

100 

92 

500 A 



1309 



Moly/Copper 

Layers 

5 X 10"^ 

115 

(809/500) 

100 

99 


* 


Mean values based on five separate tests rather than the two tests used 
for each moly/copper configuration tested. 
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total erosion depth (and in parenthesis erosion depths for the molybdenum 
or stainless steel [first] / copper [second]) measured after each test is 
recorded in column four. The durations of the tests producing these erosion 
depths is given in column five. Column six contains the calculated time 
required to erode through the depths of the two materials given in column 
four assuming they have the same erosion rates as the pure materials (from 
Figures 17 and 19). Typically the calculated erosion time agrees with the 
measured test duration to within about 10%. This suggests that multilayered 
molybdenum/copper test patches could be used in argon thrusters to determine 
erosion rates to within an accuracy of about t 10% also. 

Table VI contains similar multilayered erosion test patch data 
obtained using 105 eV mercury ions as the bombarding species. As in Table V 
the fifth and sixth columns contain the time required to conduct the test 
and the time that should have been required to erode through the measured 
depth if erosion rates for the pure base materials at the specified energy 
and nitrogen partial pressure (from Figures l8 and 20) had been used. In 
these tests the calculated erosion times are observed to be 10 to 25% below 
the measured erosion times. This suggests that a phenomenon such as coning 
may be coming into play when mercury sputtering occurs and that it may be 
reducing the erosion rates below those expected from measurements obtained 
on pure materials. Based on these results one would expect 10 to 25% 
errors in erosion rates obtained using multilayer erosion badges in mercury 


ion thrusters. 
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TABLE VI 

Sputter Erosion Characteristics 
of Multilayered Test Samples in a Mercury Ion Beam 


Test 
Samp! e 

Configuration 

Nitrogen 

Pressure 

(torr) 

Mercury 
Ion Energy 
(eV) 

Erosion 

Depth 

(A) 

Test 

Duration 

(min) 

Calculated 
Erosion Time 
(min) 

500A 

Moly/Copper 

layers 

<1 X 10"® 

105 

21 50 

(1150/1000) 

180 

163 

500A 

Moly/Copper 

layers 

5 X 10"^ 

105 

1650 

(1000/650) 

180 

142 

500A 

Moly/Copper 

layers 

-4 

1X10^ 

105 

1760 

(1000/760) 

240 

182 

500A 

Moly/Copper 

layers 

2 X 10"^ 

105 

1680 

(1000/680) 

240 

202 
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SERT II THRUSTER TESTS 

22 

Two 15 cm dia., mercury, SERT II thrusters, which were placed in 

earth orbit February 3, 1970, have been the subject of various tests since 

that time. Some of the tests have been formulated since the launch date 

as a result of the changing status of the spacecraft and thruster capabili- 

ties identified as desirable in the intervening time. It has been con- 
sidered desirable to know some characteristics of the thruster plasma and 
discharge chamber that are not measured on the spacecraft. Ground based 
tests have therefore been performed under this grant to provide data in 
support of these space tests. 

The SERT II thrusters have both been operated with their grids shorted 
together and without any high voltage applied to them. The propellant flow 
rate required to effect discharge operation at the associated currents and 
voltages were however not known. Flow rate data were measured in ground 
tests at discharge conditions similar to those measured in space and the 
results are given in Table VII. 


Table VII 

SERT II Thruster Flow Rate Data (High Voltage Off) 


Dischg. 

Dischg . 

Keeper 

Main Vaporizer 

Current 

Vol tage 

Vol tage 

Flow rate 

(A) 

(V) 

(V) 

(mA eq.) 

2.17 

39.6 

9.4 

129 

2.11 

36.6 

8.9 

153 

2.00 

34.7 

9.1 

190 


Mercury flow rates given in Table VII were measured by timing the fall 
of a mercury column in a 0.9 mm bore glass tube over a period of several 
minutes. They are considered accurate to within t 5 mA. The data of 
Table VII were obtained with the cathode heater operating at the low power 
level associated with SERT II thruster operation. Operation at the discharge 
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conditions of Table VII with this power increased to near its maximuin 
value during another test resulted in flow rates that were a few milliamps 
lower. 

The current drawn to the SERT II neutralizer keeper with the neutralizer 
off was measured as a function of keeper voltage over a range of flow rates 
first with the high voltage off and then with it on. The results of these 
measurements are given in Figure 21 along with corresponding mass flow rates 
(m) , discharge voltages (V^), and when the high voltage is on, beam currents 
(Iq). Keeper voltage is measured relative to vacuum facility ground. 

D 

With the high voltage on the collected current is seen to be smaller and 
with the high voltage off the collected current is of the order of the 65 mA 
measured on-board the spacecraft when its thrusters are operated with a 
235 volt keeper bias. The collected current is observed to increase as 
propellant flow rate is reduced with no high voltage. This behavior is 
probably a result of the increase in primary electron energy associated 
with the increase in discharge voltage brought on in turn by the flow rate 
reduction. This increase in keeper current occurs then because a high 
primary electron energy facilitates electron escape from the discharge 
chamber through the grids. 

Plasma property data were collected in the beam of a SERT II thruster 
operating at conditions close to those of the thrusters operating in space 
with the high voltage on. Figure 22 shows the variation in plasma potential, 
electron density and electron temperature as a function of radius measured 
from the thruster centerline at axial locations 5 cm, 15 cm and 75 cm down- 
stream of the accelerator grid. These data were obtained using a spherical, 
platinum Langmuir probe 2.2 mm in diameter. The thruster was operating 
under throttled conditions at a beam current of 85 mA, a neutralizer 
keeper current and voltage of 0.2 A and 28 v respectively and main, cathode 
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Fig. 21. SERT II Neutralizer Keeper Current Collection 
Characteristic - No Neutralizer Flow 
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PLASMA POTENTIAL Iv) 





Fig. 22. SERT II Beam Plasma Properties 
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and neutralizer flow rates of 11, 120 and 25 mA respectively. The power 
supplies were connected in the manner indicated in Figure 23 during the 
collection of these data. The zener diode clamping circuit acts as an 
open circuit unless the voltage across it rises above the 'v200v set point 
where it conducts. 

For a neutralizer bias voltage of zero, Vg was -9v (neutralizer cathode 
9 volts below tank ground). As the bias voltage was varied over a range of 
several volts, Vg changed such that the neutralizer cathode remained at this 
same -9v potential and the plasma properties in the beam remained unaltered. 
The plasma properties in the ion beam were also the same as those shown in 
Figure 22 when the zener diode clamping circuit was removed and replaced by 
a power supply used to hold the neutralizer bias voltage at the -9v bias 
condition where the neutralizer emission balanced the ion beam current. 
Langmuir probe trace noise was however less significant when the zener diode 
clamping circuit was in use and the neutralizer emission current was auto- 
matically forced to be equal to the beam current. 
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ELECTROTHERMAL RAMJET 

The increased use of space for the benefit of society has pointed the 
need for low cost launch systems. One class of these systems, that could 
be used to launch such materials as those needed for orbital construction 
projects or radioaactive or toxic chemical wastes for extra terrestrial 
disposal, would involve high acceleration levels. Included in this class 
of launch systems are those that involve such devices as the rail gun, mass 
driver and other electromagnetic accelerators. Because of the high accel- 
eration levels at which these devices could operate, the associated launch 
tracks would be relatively short and the terminal launch velocity could be 
achieved while the payload was still close to the earth's surface. Achieving 
terminal velocity in a short distance is considered desirable because the 
launch could be effected in a tube of reasonable size which would be sealed 
during the launch sequence if it was necessary to abort the launch. This 
could facilitate containment if the launch of toxic or radioactive materials 
were involved. 

Another propulsion concept which may be useful in these high accelera- 

23 

tion level missions is the electrothermal rocket concept. In these de- 
vices on-board propellant is heated electrically and then discharged 
through a nozzle to produce thrust. Because the heating is accomplished 
electrically, hydrogen propellant with its attendant high specific impulse 
capability can be selected. There are several types of electrothermal 
rockets. In the resistojet for example the propellant is heated by passing 
it over current carrying, resistive elements and specific impluses as high 
as 850 sec (8500 m/sec exhaust velocity) have been achieved. Higher 
specific impulses can be realized if the material limitations imposed by 
the heating elements are eliminated and joule heating is accomplished by 
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passing current directly through the propellant. These electrothermal 
rockets, known as arcjets, have demonstrated specific impulses as high as 
1500 sec (15,000 m/sec exhaust velocity). 

A basic electrothermal rocket vehicle which might be used in this 
application is shown schematically in Figure 24a. This rocket would carry 
on-board propellant which would be metered into the heating chamber where 
its temperature could be raised by energy either beamed there electro- 
magnetically or carried there from rails in the launch tube in the form 
of electrical current. 

For the mission under consideration here however where the launch is 
accomplished close to the earth's surface and in a relatively short 
distance it is also practical to distribute the propellant along the 
launch tube. When this is done the propellant can be collected by the 
launch vehicle and directed into the heating chamber as it moves along. 

The class of engine that collects its own propellant as it moves through 
its environment is the jet engine and the simplest of these is the ramjet. 
Such an engine is considered in this study because of the associated 
potential for saving energy, simplifying the design and reducing payload 
delivery costs. Physically the ramjet would operate in the manner sug- 
gested by Figure 24b by simply drawing propellant in from the tube through 
an intake diffuser, heating it electrically and then expelling it through 
a nozzle. Again energy could be transported into the heating chamber 
either el ectromagnetical ly or in the form of electrical current. No com- 
pressor or other rotating equipment would be needed although an initial 
velocity would have to be imparted to the vehicle by a booster rocket, an 
explosive charge or a low terminal velocity electromagnetic driver. 

The purpose of this study is to determine for the high acceleration 
launch mission under consideration the degree to which the performance 
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of the ideal electromagnetic accelerator might be approached by simpler, 
presumably less costly and more reliable hydrogen rocket and ramjet engines. 
The intent of the work described here is to make a preliminary performance 
comparison of several of these options for a particular mission. The 
mission selected was a high acceleration earth escape in which sufficient 
velocity would be imparted to the payload at the earth's surface so that 
it would be at earth escape velocity when it reached an altitude where 
atmospheric drag was negligible. The following propulsion schemes are 
compared in this work: 

Thruster Exhaust Velocity 

Chemical Rocket 4,500 m/sec 

Resistojet Rocket 8,500 m/sec 

Arcjet Rocket 15,000 m/sec 

Resisto-ramjet 8,500 m/sec 

Arc-ramjet 15,000 m/sec 

Electromagnetic gun 

Constraints at End of Launch Phase 

The velocity required at the surface of the earth by a payload in- 
tended to escape the earth and the maximum acceleration that such a payload 
can endure are required as inputs to the problem under consideration. These 
quantities can be defined if one considers the expulsion of the payload to 
consist of two phases, namely the launch phase illustrated in Figure 25 and 
the subsequent phase of vehicle passage through the atmosphere. A vertical 
launch trajectory is specified and initially the vehicle is assumed to be 
at rest and to have a mass m^. . When the thruster is started it generates 
a constant thrust F which produces an initial acceleration a^. . The launch 
phase is assumed to continue until time t^ when a terminal launch velocity 
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is reached. At this point the vehicle mass is m^, the acceleration is 
and the axial location measured from the initial location is x^. 
Thrusting stops and the phase of vehicle passage through the atmosphere 
begins at this point. During this second phase the vehicle velocity de- 
creases as a result of atmospheric drag. 

The maximum acceleration (a^) that the vehicle can withstand is de- 
termined by its mechanical strength. If one assumes the basic vehicle is 
essentially cylindrical and has a diameter "d" and length "1" then the 
compressive stress (S) at the base of the vehicle is equal to the thrust 

force (F) divided by its cross sectional area. 

4m. a. g 

4F ' ^ , 

S = — = (1 

ud^ Trd^ 


In this equation m^ is the final vehicle mass, is the final 
(maximum) accel eration (in G's) and g is the acceleration due to gravity. 
The ballistic coefficient of the vehicle (B), which pertains to movement 
through the atmosphere after launch, is given by 


B 


C 



(17) 


In this equation is the vehicle drag coefficient and it should have a 
value near 0.1 for the velocities under consideration here. Combining 
Eqs. 16 and 17 and solving for the maximum acceleration one obtains 


0f 


BS 

S3 * 


(18) 
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Using Figure 26 which was obtained from NASA one determines that a 

-5 2 

reasonable value for the ballistic coefficient is B = 2 X 10 m /kg. At 
this value the ratio of atmospheric entrance velocity to the earth escape 
velocity (v =11,200 m/sec) is just under 1.2 and the terminal launch 
velocity (v^) woul d therefore have to be about 1 3,000 m/sec. If one specifies an 
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BALLISTIC COEFFICIENT - B ( m*/ kg ) 


Fig. 26. Velocity Loss Through Atmospheric Drag 
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9 2 

allowable stress 1n the vehicle of 1.4 X 10 nt/m ('x-200,000 psi) then 
the maximum allowable acceleration (Eq. 18) becomes 

a. = ^ ^ ° ^ ^ = 30,000 G's 

^ (0.1) (9.8) 

This acceleration should not be exceeded if compressive stresses are to be 

9 2 

held below 1.4 X 10 nt/m during the launch to a velocity of 13 km/sec. 
This velocity will be sufficient to enable escape from the earth when the 
vehicle diameter is given by Eq. 16. 

Results 

A mathematical model developed from a dynamical force balance and 
energy conservation analysis for each of the thruster cases considered is 
presented in Appendix B. In each case constant thrust and, for the rockets 
and ramjets, constant exhaust velocity has been assumed. The results were 
obtained for a particular case where a 10 kg final mass was launched 
vertically on an earth escape trajectory. The 10 kg final mass was se- 
lected rather arbitrarily but the basic trends and comparisons made here 
are not affected by this selection. 

Figure 27 shows the velocity profiles for the six thruster options 
being compared. The effect of increasing rocket exhaust velocity from 
4,500 to 15,000 m/sec is observed to be a three-to-fourfold reduction in 
the launch time and a two-to-threefold reduction in the required launch 
track length (x^). This points out the great benefit of using the high 
exhaust velocity electrothermal rocket rather than the conventional 
chemical rocket motor which is characterized by the lower exhaust velocity. 
By using a ramjet rather than a rocket further reductions in launch times 
and track lengths are realized at a given exhaust velocity. This occurs 
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CONSTANT THRUST 
CONSTANT EXHAUST VELOCITY 
Of = 30000 G's 



Fig. 27. Velocity Profile Comparison 
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because the on-board propellant need not be accelerated in the case of 
the ramjet. At an exhaust velocity of 15,000 m/sec the ramjet performs 
essentially the same as the electromagnetic accelerator in that it requires 
the same track length (285 m) and launch time (44 msec). The 8,500 m/sec 
ramjet is seen to exhibit performance comparable to the 15,000 m/sec 
rocket. 

Figure 28 shows the time profile of the vehicle acceleration for the 
various thrusters. The acceleration is limited to 30,000 G's and this 
implies considerably lower initial accelerations for the conventional 
rockets in which the vehicle mass must decrease with time as on-board 
propellant is discharged. The shorter launch times and track lengths 
given in Figure 27 are seen to occur for the cases where the acceleration 
remains closer to the limiting value throughout the launch. 

A discussion of the variation in acceleration for the 8.5 km/sec 
ramjet will facilitate an understanding of some of the physical character- 
istics of this thruster. During the first ten milliseconds or so the 
acceleration of the 8.5 km/sec ramjet increases because on-board mass is 
being expelled to accelerate the vehicle to the point where the ram pressure 
is sufficient to effect ramjet operation. This decrease in mass causes the 
acceleration level to increase with time. During the interval from 10 to 
50 milliseconds, ramjet operation occurs and the vehicle acceleration is 
constant. Near the 50 millisecond point the vehicle velocity gets suf- 
ficiently close to the thruster exhaust velocity so that ramjet operation 

'fc 

is no longer possible. This means on-board propellant must again be used and 

* Application of the momentum balance equation demonstrates that the thrust 
(F) goes to zero as the vehicle velocity (v) approaches the exhaust 
velocity (U), unless the propellant flow rate (m) becomes infinitely 
large (F = m [U-v]) . 
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CONSTANT THRUST 
CONSTANT EXHAUST VELOCITY 
<Tf = 30000 G's 



Fig. 28. Acceleration Profile Comparison 
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acceleration again increases with time. This use of on-board propellant 
continues until the desired final velocity (13 km/sec) is realized. For 
the 15 km/sec ramjet the final phase wherein on-board propellant is used 
is not required because the final velocity is significantly less than the 
exhaust velocity. 

Figure 29 shows the variation in vehicle mass with time for each of 
the thrusters. These curves simply reflect the variations in acceleration 
observed in Figure 28. It is noteworthy however that the 4,500 m/sec 
conventional rocket requires an initial mass of 180 kg to deliver a final 
mass of 10 kg at 13 km/sec while the 15,000 m/sec ramjet can do it with 
11.5 kg and the electromagnetic thruster requires only the desired 10 kg 
as an initial mass. 

The power profiles for each operating condition are compared in 
Figure 30. The total energy (E) required to accelerate the payload to the 
specified final velocity for each option is obtained by integrating under 
the power curves and it is also shown in this figure. The variation in 
total energy required for the launch is seen to be relatively modest between 
the various modes, but the ramjet and electromagnetic designs are seen to 
require substantially greater peak power levels. Figure 30 shows that the 
peak power levels of the high exhaust velocity ramjet and the electro- 
magnetic gun are comparable. The electromagnetic gun is observed to 
require the least energy of any option, but it is important to note that 
this is a highly idealized model. The energy stored in the magnetic fields 
of this device, when the vehicle leaves the launch track, as well as joule 
heating losses will be substantial and they have not been included here. 
Their inclusion could well make the energy required for the electromagnetic 
gun more comparable to values obtained with the rockets and ramjets. 
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Fig. 29. Mass Profile Comparison 
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Fig. 30. Power Profile Comparison 
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It is important to point out that all of the modelling used here is 
highly idealized. For example losses due to aerodynamic and frictional 
drag have been neglected for the launch phase, as have propellant pumping 
powers and losses associated with energy transport through conductors or 
wave guides. These could change absolute numbers significantly although 
the comparative results are considered correct. 

Concl usions 

By using hydrogen propellant in conjunction with an electrothermal 
thruster one can increase the specific impulse by a factor of two or three 
over that obtainable with chemical systems. For a high acceleration, escape 
trajectory launch system the use of such an electric rocket translates 
into a corresponding theoretical reduction in launch times and track 
lengths of about the same factor of two or three. Further improvements 
are realized if a ramjet is used in place of a rocket. For the particular 
mission considered here an electrothermal ramjet is predicted to exhibit 
essentially the same performance as such electromagnetic guns as the rail 
gun or mass driver. Since the electrothermal ramjet has an inherent 
simplicity that makes it potentially more reliable and less costly than 
the electromagnetic devices, it appears that it could be a viable alternative. 
More detailed analysis of the electrothermal ramjet is recommended. 
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ION BEAMLET VECTORING BY GRID TRANSLATION 
John Homa 

The study of grid translation effects on ion beam optics has two 
important applications in the design of electric thrusters. For a dished 
grid pair, hole misalignments are used to redirect the emerging beamlets 
in order to minimize the off-axis thrust. A study that would yield cor- 
relations or a model of this phenomenon would facilitate design of these 

25 

grids. In addition, for some applications a vectorable beam is desired, 

and grid translation is an easily accomplished method of producing a 

moderate thrust vectoring capability. 

26 

Work done by Aston has shown the effects of such parameters as grid 
separation, grid hole diameter and net-to-total accelerating voltage ratio 
on optical performance of a two-grid system. Optical performance in this 
context is defined as the beamlet divergence variation with beam current 
as well as the maximum perveance that can be obtained. The present study 
of grid translation is intended to examine the effects of these parameters 
on vectored beams. Of particular importance is the degradation of optical 
performance as a function of relative grid displacement expressed in terms 
of beamlet divergence and maximum perveance capability of vectored beamlets 
compared to nonvectored ones. 

For a simple two-grid system the beam deflection angle can be esti- 

27 28 

mated using linear optics theory ’ by 

e = -e/(4Jig) . (19) 

where e is the relative accelerator grid displacement and 5,^ is the grid 
separation distance. The negative sign indicates that the beam deflection 
is directed opposite to the direction of relative accelerator grid 
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displacement. The parameters of Eq. 19 are defined in Fig. 31. This 
figure also illustrates pertinent grid dimensions including screen and 
accelerator grid aperture diameters (d and d ), grid thicknesses 

S d 

(t and t ), the grid separation distance [i ) and the grid offset angle 6. 
s a g 

This angle is defined in terms of the linear grid displacement e by the 
equation 

' 20 ) 


Apparatus and Procedure 

The beam vectoring study has been initiated on an 8 cm dia., mildly 

2 

divergent ion source in which the grids have been masked down to a 1 cm 
aperture area. This masking of the beam results in very uniform plasma 
properties at each aperture and hence a uniform beam. The thruster uses a 
tungsten filament cathode and neutralizer and operates with argon as the 
propellant. 

The grids are fabricated from graphite sheets with a hexagonal nineteen 

hole array. Spacing of the grids is provided by mica insulating sheets. The 

screen grid can be translated along one axis during thruster operation by a 

mechanical system coupled to a micrometer assembly as illustrated in 

Fig. 32 > while the accelerator grid is held stationary. The screen grid was 

selected as the translating grid, because such an arrangement facilitated 

grid spacing and accelerator grid alterations. 

The beam emerging from the grids is probed by a collection of molybdenum, 

2 

flat plate probes each 1 cm in area, in the manner suggested in Fig. 32. 

Twenty of these probes are arranged in a horizontal line 14 cm downstream 
from the grids. In addition, a line of probes extends vertically downward 
from each end of the horizontal probe rake. By using these three probe rakes 
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it is possible to intercept the full beam profile over the complete range 
of grid translation. A 70% open area screen is placed in front of the 
probes, and is biased 30 v negative relative to ground in order to reflect 
neutralizer electrons. The probes themselves are biased 18 v positive 
relative to ground to reflect low energy, charge-exchange ions produced 
downstream of the grids. The rake assembly is designed so it can be 
rotated in the manner suggested in Fig. 32 to facilitate rake axis align- 


ment with the beam axis in the plane perpendicular to the direction of grid 
translation. Data was initially collected using only the horizontal probe 
rake, however beam profiles for vectored beams showed significant currents 
being collected on the end probes. The vertical probe rakes were then 
added to insure that the entire beam profile was obtained. 

Standard thruster operating conditions for preliminary checkout have 


been a total accelerating voltage of = 600 v with a net accelerating 

voltage set to obtain R = values of 0.5, 0.7, and 0.9. The 


discharge voltage has been 40 v and operating pressures in the range of 

_5 

2-4X10 torr have been used. Once a beam current density profile has 
been measured across the beam diameter it is analyzed using a computer 
routine. This routine first uses the current density data at the probe 
rake locations to compute the current density on a circular arc located on 
a 17 cm from the grid center. From this profile the point of maximum 


current density is identified. The deflection angle is then computed as 
the angle between the perpendicular to the grids and the line from the 


screen grid center through the point of maximum current density. Next the 


two segments of the current density profile on either side of the maximum 
are averaged and the total beam current is estimated by integrating this 
average profile for an assumed axi-symmetric distribution about the line 
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of maximum current density. Inherent in this portion of the analysis is 
the assumption that the ions are drawn from a point source; this assumption 
does not introduce a significant error because the beam diameter is small 
compared to the distance to the probe rake. The beam current obtained in 
the analysis by integrating the current density is compared to the 
measured beam current to insure that they agree to within 5-10%. The 
divergence half angle (a in Fig. 31) is computed as the half angle of the 
cone enclosing 95% of the integrated beam current determined in the afore- 
mentioned analysis. This definition of half angle is the same as the one 

26 

used in previous work. The 95% criterion was selected as a compromise 
value that yielded a sensitive indicator of ion beam deflection. A thrust 
coefficient defined as the fraction of the integrated current that is 
directed along the deflected beam axis is also calculated. 

Resul ts 

The data that have been collected so far during this study are con- 
sidered preliminary, but they do seem to be self-consistent and in agreement 
with results obtained by others. Figure 33 shows typical results wherein 
the beam deflection angle (g) is seen to vary linearly with the grid 
offset angle (6) over the range 1 15° (+ 10° deflection angle). Beyond 
this point, the accelerator grid impingement current increases rapidly, due 
to direct interception of the beam, and the beam deflection becomes non- 
1 inear. 

The divergence angle and thrust coefficient are shown in Fig. 34, as 
a function of grid displacement. The graph shows a general increase in 
divergence angle which reaches a peak near the point at which accelerator 
impingement increases. The decrease in divergence angle with further grid 
displacement is probably due to interception by the accelerator grid of 




Fig. 33. Beamlet Deflection Characteristic 
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the beam edge resulting in a less divergent beam. It is noteworthy that 
the divergence angle is not a strong function of grid offset angle for 
the particular grid and operating condition being investigated here. 

Conclusion 

The beam vectoring apparatus functions properly and consistent results 
have been obtained. Further work is needed to investigate the effects of 
variations in accelerator and screen grid thicknesses and hole diameters, 
grid separation total accelerating voltage, and net-to-total accelerating 
voltage ratio. 
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MULTIPLE CATHODE STUDY 
John R. Brophy 

Introduction 

One method of increasing the thrust density from electrostatic ion 
engines is to increase the ion beam current density extracted by the 
accelerator system. This in turn requires an increase in the plasma 
density from which the ions are extracted. To increase this plasma density 
may require that the thruster be operated at discharge currents of several 
tens of amperes. Cathodes operating at these current levels may also be 
required in large diameter thrusters. Two methods of generating these high 
discharge current levels have been suggested. The first requires the 
development of a high power cathode capable of emitting tens of amperes 
while still having a lifetime comparable to those of existing cathodes. The 
second method would involve parallel operation of several conventional hollow 
cathodes. The number of cathodes would in this case be selected so the 
emission current requirement for each individual cathode would be within 
the capability of presently available designs. The problem with this 
latter solution is the tendency for cathodes to share the emission load 
unequally during operation or for some of the cathodes to fail to start. 

The experiment described in this section was designed to determine how 
this potential problem might be overcome. 

Apparatus 

29 

A 15 cm diameter multipole thruster was modified to operate with 
two hollow cathodes supplying electrons to the main discharge chamber in 
the manner suggested by Fig. 35. The cathode discharge region as shown 
here is seen to be a torus of rectangular cross section that is connected 
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to the main discharge chamber by an annular baffle aperture on its down- 
stream side. The cathodes are located diametrically opposite each other at 
the upstream side of the cathode discharge region. Each cathode is elec- 
trically isolated from the other one and from the thruster body as suggested 
by the electrical schematic shown in Fig. 36. In order to achieve complete 
isolation of the two cathodes it was necessary to use separate keeper, 
vaporizer and tip heater power supplies as well as separate propellent 
feed systems for each cathode. Electrical isolation of the cathodes was 
required to facilitate measurement and subsequent control of the two 
emission currents. The separate propellent feed systems were required not 
only for the electrical isolation of the cathodes but also to facilitate 
independent measurement of the flow rates to each cathode. 

Discussion 

The first problem encountered in the operation of the multi-cathode 
thruster was associated with the ignition of both hollow cathode/keeper 
discharges. On several occasions it was observed that the initiation of 
one of the cathode/keeper discharges, while the second discharge failed 
to start, resulted in coupling of the plasma from the operating cathode 
to the keeper of the nonemitting cathode. This coupling effectively 
prevented the ignition of the second cathode. To assure that both cathodes 
could be started it was found that they had to be isolated from each other, 
either with a physical or a magnetic barrier, that prevented this coupling. 
In the thruster configuration tested here it was found that increasing 
the magnetic baffle current would create a leakage magnetic field of 
sufficient magnitude to isolate the cathodes. Isolating the cathodes with 
a magnetic field in this manner generally allowed both cathodes to be 
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started without great difficulty. 

After both cathodes had been started and the main discharge had been 
initiated the problem became one of control of the relative emission 
currents from the two cathodes so they were nearly equal. It was found 
that this control could be effected by placing variable resistors in the 
leads to each cathode in the manner suggested by Fig. 37. The resistance 
in each lead could then be adjusted until the emission currents were equal. 
Figure 38 indicates the magnitude of the resistances required to equalize 
the emission currents at several discharge current levels. Emission 
currents shown in Figure 38 have been normalized with the discharge current. 
For the cathodes used in this study it was found that the one designated 
as the left cathode would, under normal operating conditions, tend to emit 
most of the required current. For this reason the resistance in the lead 
to the right cathode was set to zero and the resistance in the lead for the 
left cathode was varied from zero to lOfi to generate the data shown in 
Fig. 38. Two interesting things should be noted from this Figure. First, 
the resistance required is independent of the magnitude of the discharge 
current, and second, the total cathode emission current is less than the 
current collected at the anode (i.e. the discharge current). In fact the 
emission current in these cases accounted for just under 70% of the dis- 
charge current. The remaining 30% of the discharge current is due to the 
electrons left behind by ions which either recombine on the discharge 
chamber walls or are extracted into the ion beam. The details of the 
main discharge current balance are discussed in Ref. 2. 

One disadvantage of using resistances to control the relative emission 
between cathodes operating in parallel is that power is wasted in the 
resistors. For this reason control of the relative emission currents was 
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FOR EMISSION CONTROL 


Fig. 37. Location of Tuning Resistors for Control of 
Relative Cathode Emission 






NORMALIZED CATHODE EMISSION, U/I 
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Fig. 38. Cathode Emission Control Using Tuning Resistors 
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attempted through adjustment of the cathode tip heater powers and neutral 
flow rates only. It was found that the response time of the relative 
emission currents to manual adjustments of these parameters was sufficiently 
rapid to allow control to be achieved in this manner. Typical test results 
showing the emissions from two parallel cathodes together with variations 
in cathode heater currents and cathode propellant flow rates required to 
achieve control are shown in Fig. 39. Initially the left cathode had an 
emission current of 3.8 amps and the right cathode was emitting 2.0 amps as 
shown in Fig. 39a. By turning the left cathode tip heater off and adjusting 
the right cathode heater current and neutral flow rates during the test in 
the manner shown in Figs. 39b and 39c, the emission currents were maintained 
essentially equal from the time 30 minutes into the test to its conclusion. 
The frequent small adjustments of the heater current required to maintain 
control can be seen in Fig. 39b. The less frequent adjustments required 
in the neutral flow rates can be seen in Fig. 39c. 

Concluding Remarks 

It is possible to operate an electron bombardment ion thruster with 
multiple cathodes supplying the emission current to the main discharge 
chamber. Starting more than one cathode can be difficult if the cathodes 
are not separated either physically or with a magnetic field. The most 
effective physical separation would probably be for each cathode to be in 
its own separate cathode discharge chamber. Control of the emission current 
split between the cathodes can be achieved through the use of variable 
resistors in the leads to the cathodes and by controlling the tip heater 
powers and neutral flow rates to these cathodes. 



RIGHT CATHODE 

NEUTRAL FLOW HEATER CURRENT EMISSION CURRENT 

RATE (m A eq.) (amps) (amps) 
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o RIGHT CATHODE EMISSION 
□ LEFT CATHODE EMISSION 
o TOTAL EMISSION 




Fig. 39. Control Parameter Profiles 



107 


MODIFIED BOHM VELOCITY FOR PLASMA CONTAINING 
PRIMARY ELECTRONS AND DOUBLY CHARGED IONS 
John R. Brophy 

The Bohm criterion is useful in ion thruster work because it provides 

a model describing the rate of ion migration toward discharge chamber 

30 

boundaries. The initial model developed by Bohm assumed only Maxwellian 

31 

electrons were present in the plasma. Masek modified this to account for 
the density of primary electrons. The purpose of this work is to investi- 
gate the effects of both the energy of the primary electrons and of 
doubly-charged ions that are generally present in ion thruster discharge 
chamber plasmas. Consideration of these species is becoming increasingly 
important as thruster size and power density increase and induce a corres- 
ponding, increase in the densities of doubly-charged ions and primary 
el ectrons . 

The one-dimensional Poisson equation that applies in the sheath region is, 

0= [n^(x) + 2n^^(x) - n^(x) - n (x)] (21) 

where (|), the potential in the sheath, is taken as zero at the plasma sheath 

boundary, e is the electronic charge and is the permeability of free space. 

The position variable (x) in this equation is taken to be zero at the sheath 

edge and n., n. ,» n and n are the singly-charged ion, doubly-charged ion, 

+ ++ m p 

Maxwellian electron and primary electron densities respectively. The 
Maxwellian electron density within the sheath is determined by the baro- 
metric equation in the form, 

ei>/kT 

nj^(x) = r\^{6) e (22) 

where n^(o) is the electron density at the plasma sheath interface, T is 
the temperature of the Maxwellian population and k is the Boltzmann 
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constant. For the primary electrons a corresponding equation can be 
developed by considering the continuity of primary electrons through the 
sheath; 

np(x)Vp(x) = np(o)Vp{o)(l + (23) 

C [J 

In this equation v (x) and v (o) are the velocities of the primary electrons 

r r 

in the sheath and at the sheath edge respectively. The term multiplying 
the right hand side of this equation accounts for electron reflection in- 
duce by the electric field in the sheath. Conservation of energy for the 
primary electrons in the sheath is expressed through the equation 


^mgV2(x) - e (|)(x) =|-mgV2(o) (24) 

where m is the electron mass. Combining Eqs. 23 and 24 one obtains the 
0 

primary electron density variation through the sheath in the form. 


rip(x) = np(o) [1 



(25) 


where n^Co) is the primary electron density at x = 0 and e 
the primary electrons in the bulk plasma, is given by. 


P 


the energy of 


£p=-^Vp^°^ • 

The singly and doubly-charged ion densities in the sheath are also found by 
combining their respective continuity equations. 


n^(x) v^(x) = n^(o) v^(o) for singly-charged ions 

n^^(x) v^^(x) = n^^(o) v^^(o) for doubly-charged ions 

with the corresponding energy conservation equations, 

m^v2(x) =^m^v2(o) - e (|)(x) for singly-charged ions 

^m.v2^(x) = |-m^.v2^(o) - 2e (|)(x) 


for doubly-charged ions 


(28) 
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to obtain. 




n++(x) = n^^(o) 


for singly-charged ions 
for doubly-charged ions 


(29) 


In these equations n^(o) and n^^(o) are the singly and doubly-charged ion 
densities at x * o respectively and and are the corresponding 
energies at x = o. Substituting Eqs. 22, 25 and 29 into Eq. 21 yields. 


d2 


dx 


e_ 

e. 


n^(o) 


- e(j) 


+ 2n^^(o) 


'++ 


'++ 


2eij) 




(30) 


Multiplying Eq. 30 by d<()/dx, integrating and solving for the constant of 
integration satisfying the boundary condition d(|)/dx = 0 at (|) = 0 yields. 




4e 


n^.(o) 




^ 'V<!l 

e ^++' 


r 


M - 1 


■++ 


kin (o) 
m ' 

2e 




(31) 


But the quantity (d({)/dx)2 must always be greater than or equal to zero 
because an imaginary potential gradient is not physical ly meaningful . 
Application of this constraint yields 



This equation can be simplified by expanding the various terms using the 
following power series: 



no 




% 


0 _|±, 


(T . M) =1.^1 (£i 


(33) 


'++ 


^++ ^ 'S+ 


(1 + M.) / = 1 + M+ 3 

"p ® "p 

These expressions are now substituted into Eq. 32 and it is recognized that 


^ «1 , M «i , ^ ^ ei 

^++ 


(34) 


so one need only retain the terms in Equation 33 that are of quadratic or 
lower powers. The result is 

^2\ 

\ , . / f \ ■m"'' / 

" ?RT 

^ (35) 




2e 




Recalling that charge neutrality exists at the plasma sheath boundary the 
equation 

2n++(o) + n^(o) = n^(o) + np(o) (36) 

may be substituted into Eq. 35 to yield. 


"+(o) (o) 2n (o) n (o) 

_ 11 + -^ + _2 ;o . 

S S+ ^p 

Noting that the doubly-charged ion energy at the sheath edge should be 
twice the single ion energy, i.e. 


(37) 


£ ” 2e 

^++ 


( 38 ) 



m 


Eq. 37 becomes. 


-n+(o) 2n^^(o) 2n (o) n (o) 


m' 

kT 


> 0 . 


+ + ■■ -P 

Solving Eq. 39 for the singly-charged ion energy at the plasma sheath 
boundary yields, 

n^(o) t 2n^^(o) 

' 2n (o) n (o)\ ’ 

kT e. 


(39) 


(40) 


Making use of the relationship 


1 2 
2 


(41) 


Eq. 40 can be written in terms of the ion velocity at the sheath edge as. 


^T 

‘ (n^(o 

) + 2n^^(o))/ 


m. 

[ 1 + 

rkT \ 




(42) 


If Eq. 36 is substituted into Eq. 42 the result is. 



(43) 


This equation is the one sought initially namely the one reflecting the 
effects of primary electron energy and doubly-charged ions. For the usual 
case where 


KT 


S 1^ 


«1 


this simplifies to 

V , 


\/k7' 

n (o) 

1 r ^ . 

■V”i 

n (o) 
m 


(44) 


which is the result given by Masek 
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It is interesting to note that the 


doubly-charged ion density does not appear explicity in Eq. 43. The Bohm 



velocity for doubly-charged ions is obtained from Eq. 43 by using 


St = 2s = 2 ' “ri 


(45) 


or 


to obtain 


%+ = »t 



'2kT\/’ ^ 


,1 + 


kT 

2e 


I M^)/ 

rv m ' ' 


p m 


In the limit of a high primary-to-Maxwel 1 ian electron density ratio 
Eqs. 43 and 46 suggest ion velocities of 

V = _!± - l/_^ • 

+ ^ y ni. 

This appears to be a physically reasonable limiting condition. 


(46) 


(47) 
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APPENDIX A 

Langmuir Probe Sweep Circuitry 

The circuit diagram for the equipment designed to sweep the Langmuir 
probe over a desired voltage range is shown in Figure A1 . The sweep circuit 
employs an operational amplifier-integrator circuit to produce a voltage 
ramp. This signal is then fed into a series pass power transitor to provide 
the necessary current output. Input terminals and switching are also 
provided for biasing the probe to keep it clean. A trace is taken by first 
pressing button S^. This switches the probe from the cleaning circuit to 
the voltage sweep circuit and at the same time engages the servo for the x-y 
recorder. Then pressing button S-j activates the recorder pen and begins the 
voltage sweep. The sweep rate is controlled by adjusting resistor R^, the 
start voltage by resistor R^. and the final voltage by resistor Ry. The 
sweep rate range can also be changed by changing the value of capacitor C-| . 
The circuit will provide currents up to 500 mA for voltages swept between 
-15 to +15 V over a wide range of sweep rates. A probe cleaning circuit is 
also shown as an integral part of the design. Switches 5 and 6 provide 
manual override for the recorder pen and servo to aid in setting up the axes 
for the trace. The x-voltage follower uses an operational amplifier in a 
noninverting unity gain amplifier circuit to drive the x-channel of the 
recorder without drawing any appreciable current through the current sensing 


resistor RIO. 
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APPENDIX B 

Electrothermal Launch System Dynamics 
Nomenclature (mks units) 

F - Thrust Force 

g - Acceleration due to gravity 

• 

m - Mass flow rate of Propellant 
m - Mass of Vehicle 

P - Input Electrical Power 

t - Time 

U - Thruster exhaust velocity 

V - Vehicle velocity 

X - Position 

a - Acceleration (G's) 

Subscripts 

i - initial condition 

f - final condition 

p - pertaining to the pumping process 

If - pertaining to the end of phase one in ramjet mode 

2f - pertaining to the end of phase two in ramjet mode 

In this development it has been assumed that 

• the launch would be vertical, 

• drag and frictional forces are negligible, 

t any thruster exhaust is fully expanded isentropical ly. 
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• any engine operates at constant thrust and exhaust velocity, and 

• the launch is effected close to the earth's surface so the 
acceleration due to gravity is constant. 

The basic schematic of the launch process has been shown previously in the 
text (Figure 25 - page 78). 

Rocket Model 

The thrust of a rocket engine which receives its propellant from on-board 
tanks is given by 

F = mU. (Bl) 

Balancing this thrust against the weight and acceleration of the vehicle 
yields 

F = mga = m.ga. = gm + m^ . (B2) 

Integration of Eqn. Bl gives vehicle mass as a function of time 

m - m It 

" "’i ■ U (B3) 

and combining Eqs. B2 and B3 one obtains an expression for the vehicle 
acceleration as a function of time 



Integrating Eq. B2 one obtains the following expression for vehicle velocity 
as a function of time 

) v = -U iin(l - £ ) - g t . (B5) 

i 

A second integration of this equation gives position as a function of time. 

* 

X = (m. - Jin(l - -Oi.) + Ut - -31^ (B6) 

■# 

The input power required is given by 

P = ^ + Pp (B7) 



where is the pumping power. A preliminary analysis has shown the 
pumping power is only a few percent of the total power so it will be 
neglected here. 

Electromagnetic Accelerator Model 

For this case the vehicle mass remains constant as does the thrust 
defined by the electromagnetic mechanism used to produce it. In this case 
the force balance on the vehicle becomes 

F = m^.ga = m^.ga^. = gm^. + ^ . (B8) 

The vehicle mass is given by 

m = m^ (B9) 

and the acceleration for Eq. B8 by 

a = a^. . (BIO) 

Integration of Equation B8 gives the velocity as a function of time 

V = (^ - g)t . (BID 

i 

A second integration gives the variation in position with time 

X = (|- - g) . (B12) 

1 

The required input power in this case will be idealized to the extent 
that losses associated with magnetic energy stored at the conclusion of the 
launch phase and joule heating losses are neglected. For this ideal 
launch case the input power is given by 
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Ramjet Model 

During the first phase of operation of a ramjet it is necessary to 
increase the vehicle velocity from zero to the point required to produce 
the stagnation pressure needed for proper engine operation. This initial 
acceleration could be accomplished using an expanding explosive, an 
electromagnetic starter or an auxiliary rocket. The rocket is assumed 
here because it seems simplest although it introduces some penalty in terms 
of vehicle mass required. With the application of thrust using on-board 
propellant and a rocket motor Eqs. B1 through B7 describe the first phase 
of operation. This phase is continued until the vehicle velocity is 
sufficient to satisfy the equation 

2 P* 

^ (B14) 

P 

where P* is the required engine stagnation pressure and p is the propellant 
density. For the case under consideration here it was determined that a 
vehicle velocity about 15% of the jet velocity was sufficient to assure 
the required pressure. 

During the second phase of ramjet operation the inertial force of the 
incoming propellant must be considered and the net thrust on the vehicle 
obtained from a momentum balance becomes 

F = m (U - v) . (B15) 

During this phase the vehicle mass remains constant at a value of 

Ft 

m = m^^ = m^. g — . (B16) 

Because the mass is constant the acceleration of the vehicle also remains 
constant and is given by 
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a = 


gtn^f 


g(m^- - 


Ft 

T 


If 


(B17) 


) 


Integration of Equation B2 under the boundary condition that the 
velocity is v^^ at time t^^ and using the fact that the mass remains 
constant during this phase gives 


'' = ''if " (i 


If 


g) (t - t^^) . 


(B18) 


A second integration gives the variation in vehicle position with time 


^ ~ ^If ''if^^ ■ '*^1f^ "^2 ^ m ^^ - g) (t - 


(B19) 


The power consumption during this phase must reflect the fact that some of 
the power input to the propellant comes from the velocity of the propellant 
relative to the vehicle and this can be subtracted from the power of 

ml |2 

the exhaust stream Thus the input power becomes 

mU^ mv^ 


P = 


(B20) 


where the actual pumping power required has again been neglected. 

A third phase of ramjet operation is required when the desired 
terminal velocity is greater than the vehicle exhaust velocity. The 
reason for this is seen from Eq. B15 which shows the thrust goes to zero 
as the vehicle velocity approaches the exhaust velocity. Under the con- 
straint of operation at constant thrust and exhaust velocity being con- 
sidered here, this equation also suggests the mass flow rate would have to 
become infinite as v U. In order to limit the required mass flow rate 
to a reasonable value this third phase has been defined rather arbitrarily 
as starting when the vehicle velocity is 90% of the exhaust velocity. The 
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« 


« 


basic equations written for phase one apply here again although integration 
between different limits is required. 

The expression for the vehicle mass is obtained by integration of 
Eq. BIO. 

m = niif - (t - t2f) . (B21 ) 


The acceleration is obtained by combining Eqs. B2 and B21 to obtain 


a 


g[m-,f - -g - t 2 ^] 


(B22) 


Integration of Eq. B2 using the initial conditions for this phase gives 


V - V2^ - U Jin [1 - (t - t2^)] - g(t - t 2 ^) (B23) 

Integration of this equation yields 

m. _U P 

^ ^ ^2f ''2f^^ “ ^2f^ ^ ■ ^2f^^ ■ m^ ■ ^2f^^ 

+ t - t2f} - ■|(t - t2^)^ . (B24) 

For each mode of operation the total time required for the launch is 
determined by the equation for vehicle velocity through the specification 
of desired final velocity. 
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